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ABSTRACT 

We present a photometric survey of the optical counterparts of ultraluminous X-ray sources (ULXs) observed 
with the Hubble Space Telescope in nearby (< 5 Mpc) galaxies. Of the 33 ULXs with Hubble & Chandra 
data, 9 have no visible counterpart, placing limits on their My of ~ -4 to -9, enabling us to rule out O-type 
companions in 4 cases. The refined positions of two ULXs place them in the nucleus of their host galaxy. They 
are removed from our sample. Of the 22 remaining ULXs, 13 have one possible optical counterpart, while 
multiple are visible within the error regions of other ULXs. By calculating the number of chance coincidences, 
we estimate that 13 ±5 are the true counterparts. We attempt to constrain the nature of the companions by 
fitting the SED and My to obtain candidate spectral types. We can rule out O-type companions in 20 cases, 
while we find that one ULX (NGC 253 ULX2) excludes all OB -type companions. Fitting with X-ray irradiated 
models provides constraints on the donor star mass and radius. For 7 ULXs, we are able to impose inclination- 
dependent upper and/or lower limits on the black holes mass, if the extinction to the assumed companion star is 
not larger than the Galactic column. These are NGC 55 ULX1, NGC 253 ULX1, NGC 253 ULX2, NGC 253 
XMM6, Ho IX X-l, IC342 X-l & NGC 5204 X-l. This suggests that 10 ULXs do not have O companions, 
while none of the 18 fitted rule out B-type companions. 

Keywords: Accretion: accretion discs, Black hole physics, Stars: binaries: general, X-rays: binaries. 



1. INTRODUCTION 

In the late 1970s and early 1980s the Einstein telescope 
was used to perform studies of normal galaxies, which re- 
vealed the presence of X-ray luminous non-nuclear objects 
that were brighter than those in our own galaxy (Fabbiano 
1989). These objects were later termed ultraluminous X-ray 
sources (ULXs; e.g. Makishima et al. 2000). Although stud- 
ies of these luminous sources have continued for more than 
30 years, their nature is still unclear (e.g. Roberts 2007, Glad- 
stone 2010, Feng & Soria 2011). It has been confirmed that 
many of these sources contain accreting black holes (Kubota 
et al. 2001), but currently the masses of the compact objects 
are still unknown. Their luminosities (Lx > 10 39 erg s _1 ) 
preclude the possibility that we are observing isotropic, sub- 
Eddington, accretion onto a stellar mass black hole (sMBH; 
comparable to the Galactic sMBHs, 3 < M BH < 20M o ), 
which has led to the idea that we may be observing interme- 
diate mass black holes (IMBHs; 10 2 - 10 4 M ; Colbert & 
Mushotzky 1999). An alternative is that we may instead be 
observing massive stellar remnant black holes (MsBHs, Feng 
& Soria 201 1; defined as the end product from the death of a 
single, current generation star; Mbh ^ 1OOM ; e.g. Fryer 
& Kalogera 2001; Heger et al. 2003; Belczynski, Sadowski 
& Rasio 2004) that are either breaking or circumventing the 
Eddington limit. 
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X-ray analysis has been exhaustive, with early XMM- 
Newton observations providing both spectral and timing ev- 
idence viewed as supporting IMBHs. Analysis of their X- 
ray spectra revealed the presence of cool disc emission, in 
combination with a power-law, a combination used effec- 
tively to study Galactic sMBH systems. However, here the 
disc temperature appeared much cooler than those of sMBHs, 
suggesting IMBHs (e.g. Miller et al. 2003; Kaaret et al. 
2003; Miller, Fabian & Miller 2004). Meanwhile, X-ray tim- 
ing analysis revealed the presence of quasi-periodic oscilla- 
tions (QPOs) in M82 X-l, M82 X42.3+59 & NGC 5408 X- 
1 (Casella et al. 2008; Feng et al. 2010; Strohmayer et 
al. 2007), something observed in both Galactic sMBH and 
SMBH systems, that seems to scale with black hole mass 
(e.g. McClintock & Remillard 2006; McHardy et al. 2006; 
van der Klis 2006). Combining both the spectral and timing 
analysis of these sources indicated that ULXs were hosts to 
IMBHs (e.g. Dheeraj & Strohmayer 2012). However, recent 
timing studies of ULXs show that many of these systems ap- 
pear to have suppressed intra-observational variability (Heil 
et al. 2009), while spectral studies of higher quality data have 
indicated the presence of a break above 3 keV, a feature that 
would not be expected if we were viewing sub-Eddington ac- 
cretion onto an IMBH (e.g. Stobbart et al. 2006; Gladstone, 
Roberts & Done 2009). Instead, we have suggested that this 
combination of new X-ray spectral and timing features de- 
scribes a new super-Eddington accretion state, the so called 
ultraluminous state (Roberts 2007; Gladstone et al. 2009; 
Roberts et al. 2010). This is echoed in the re-analysis of 
NGC 5408 X-l, which shows that the spectra, power spec- 
tral density and rms spectra of this source are better matched 
to models of super-Eddington accretion than sub-Eddington 
accretion onto an IMBH (Middleton et al. 2011). In Mid- 
dleton et al. (2011a), the authors reported on the analysis of 
XMM flux-binned data from M33 X-8, and suggested a simi- 
lar two component spectral fit. The spectral evolution and tim- 
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ing properties are unlike those of the standard sub-Eddington 
accretion states, leading the authors to invoke the onset of an 
extended photosphere and a wind to explain the observed data 
(photosphere and/or outflow dominated, as predicted in mod- 
els of super-Eddington accretion; e.g Begelman et al. 2006, 
Poutanen et al. 2007). 

Although recent X-ray analysis seems to be pointing to- 
wards the presence of MsBHs, more direct evidence is re- 
quired. An attractive method is that used to confirm the first 
known Galactic black hole, Cygnus X-l. Following Mur- 
din & Webster (1971), we must first find the potential optical 
counterpart photometrically. Optical spectroscopic follow-up 
can then be performed to gain dynamical mass estimates for 
the system (Webster & Murdin 1972; Bolton 1972; Paczynski 
1974). Such techniques have been used repeatedly with suc- 
cess in our Galaxy, e.g. van Paradijs & McClintock (1995) 
and Charles & Coe (2006), and more recently on the extra- 
galactic source IC 10 X-l (Prestwich et al. 2007; Silverman 
& Filippenko 2009). If we could apply this method to ULXs, 
we could settle the debate over the mass of the black holes 
contained within these systems. 

Here, we use optical observations to identify and classify 
the optical counterparts to ULXs. The identification of unique 
counterparts is not trivial, as many of these objects reside in 
crowded stellar fields (Liu et al. 2009), which is unsurpris- 
ing given their apparent association with star-forming regions 
(Fabbiano, Zezas & Murray 2001; Lira et al. 2002; Gao et 
al 2003; Swartz, Tennant & Soria 2009). This suggests that 
we may be looking for young companions, a theory supported 
by recent results finding a prevalence for blue companions to 
these ULXs indicative of OB-type stars (e.g. Liu et al. 2004; 
Grise et al. 2005; Roberts et al. 2008), but it should be noted 
that such a blue colour may be partly due to contamination 
by reprocessed X-rays from the accretion disc or stellar sur- 
face (e.g. Copperwheat et al. 2005; 2007; Madhusudhan et al. 
2008; Patruno & Zampieri 2010; Grise et al. 2012). 

To date much of the analysis of potential counterparts has 
focused on individual source and/or their host galaxies (e.g. 
Ho II X-l, Kaaret, Ward & Zezas 2004; NGC 1313 X-2 Muc- 
ciarelli et al. 2005; Liu et al. 2007; Grise et al. 2008; Im- 
piombato et al. 2011; NGC 5408 X-l, Lang et al. 2007; 
M51 population, Terashima, Inoue &Wilson 2006; Antennae 
galaxy, Zezas et al. 2002; Cartwheel galaxy, Gao et al. 2003), 
whilst only a small number of larger surveys have taken place 
(e.g. Ptak et al 2006, Swarz, Tennant & Soria 2009; Tao et al. 
2011). 

Spectroscopic follow-up has begun for some of these 
sources, although only a small number have been published 
to date. Roberts et al. (2001) studied the counterpart of NGC 
5204 X-l, finding a blue, almost featureless spectrum, with 
similar featureless spectra found for other sources (e.g. NGC 
1313 X-2, Zampieri et al. 2004; Roberts et al. 2011; Ho IX 
X-l, Grise et al. 2011; Roberts et al. 2011; NGC 5408 X-l, 
Kaaret & Corbel 2009; Cseh et al. 2011; Grise et al. 2012). 
This suggests the light is non-stellar in origin (see Figure 1 in 
Roberts et al. 2011; Gladstone et al. in prep), indicating that 
the light may be dominated by emission from the accretion 
disc. Nevertheless, the search for dynamical mass constraints 
has continued, as a number of spectra contain the He II 4686 A 
high excitation line. This line has been associated with accre- 
tion discs in Galactic sources, and has been used successfully 
to gain such mass constraints in the past (e.g. GRO J 165 5 -40; 
Soria et al. 1998). Initial results from the optical analysis of 



multi-epoch spectra of two ULXs (NGC 1313 X-2 & Ho IX 
X-l) have detected radial velocity variations; however, they 
may not be sinusoidal (Roberts et al. 201 1; Gladstone et al. in 
prep). Studies have also been performed on the optical coun- 
terpart to ULX P13, in NGC 7793 (Motch et al. 2011). Here 
variations in the He II line are also present, but superimposed 
on a photospheric spectrum. This reveals the possible pres- 
ence of a late B-type supergiant companion of between 10 - 
20 M (Motch et al. 2011). Again, radial velocity variations 
were detected for this source, with further data and analysis 
required to confirm its period and nature. 

To date the confirmation of black hole mass has proved elu- 
sive for known counterparts. This paper seeks to find more 
potential counterparts for further study. We focus our search 
on nearby galaxies, in order to have the best chance of finding 
unique optical counterparts, while maximising the potential 
for photometric and spectroscopic follow-up of these systems. 

The paper is arranged as follows. First we outline the sam- 
ple selection, and the data reduction processes (Sections [2] & 
|3]). We then go on to combine optical and X-ray imaging data 
to identify all possible counterparts in Section ffl Section [5] 
applies multiple techniques to classify these candidate coun- 
terparts, while Section [6] presents a discussion of our results, 
implications of the analysis, and routes for further study. 

2. SAMPLE SELECTION 

We compiled a complete list of the known ULXs, draw- 
ing from a number of ULX catalogues available in the public 
domain, including Roberts & Warwick (2000), Swartz et al. 
(2004), Liu & Mirabel (2005), Liu & Bregman (2005), Ptak 
et al. (2006) and Winter, Mushotzky, & Reynolds (2006). A 
concise primary list was formed by merging duplicate identi- 
fications and removing any sources for which subsequent re- 
search has indicated that a ULX was not present (based on 
luminosity criteria, or the object later being identified as a 
non-ULX, e.g. a foreground star or background quasar). 

Many of the optical counterparts identified to date are faint 
(> 24 mag; Roberts, Levan & Goad 2008). Therefore, 
we place an additional distance constraint on our sample of 
5 Mpc, to allow for potential photometric & spectroscopic 
follow-up. At this distance, a B0 V star would have an appar- 
ent magnitude of 24.4 (My = -4.1, Zombeck 1990), so more 
distant objects would be impractical for spectroscopic studies 
with current international ground based facilities. 

We retain the ULXs residing within NGC 3034 (M82) in 
our sample as some distance estimates have indicated that this 
galaxy may be located within 5 Mpc (e.g. ~ 3.6 Mpc, Freed- 
man et al. 1994). This provides a sample of 45 nearby ULXs 
that we list in Table [T] along with their published luminosi- 
ties, distances, Galactic absorption and extinction columns 
for both the optical and X-ray bands. We also include extra- 
galactic or total Nh columns for each source, as found via 
literature search (where available). It is not clear where the 
additional absorbing material is located, as it could be gas 
clouds in the host galaxy, or associated with the ULX itself 
(e.g. photosphere/wind). 

Table [T] also indicates that the sources residing within 5 
Mpc cover the majority of the standard ULX X-ray luminos- 
ity range (Lx ~ 10 39 - ~ 10 41 erg s -1 ; only excluding the 
new hyperluminous X-ray sources e.g. Sutton et al. 2012). 
The X-ray luminosities listed within this table are taken from 
the references denoted in superscript. As the data is collated 
from published results we highlight inconsistencies in their 
calculation. Many are the observed X-ray luminosities, al- 
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though some are intrinsic/de-absorbed (identified by T). The 
luminosities listed in Table Q] include values derived from ob- 
servations using three separate X-ray telescopes (ROSAT - R, 
Chandra - C and XMM-Newton - X), with detectors that are 
sensitive to differing energy ranges. 

Liu & Bregman (2005) used the ROSAT archive in combi- 
nation with the online tool WEB PIMM ^]to extrapolate a lumi- 
nosity over the 0.3 - 8.0 keV energy range. In each case they 
used a photon index of 1.7 and Galactic Nh- Luminosities 
calculated using webpimms are marked with a 'P'. The au- 
thors Swartz et al. (2004), Humphrey et al. (2003) and Bauer 
et al. (2001) each used data from the Chandra X-ray tele- 
scope to provide luminosities over the ranges 0.5 - 8.0, 0.3 

- 7.0 and 0.5 - 10.0 keV ranges respectively. Some authors 
provide 0.3-10 keV luminosities calculated using observa- 
tions from the XMM-Newton telescope including Winter et al. 
(2006), Stobbart et al. (2006) and Feng & Kaaret (2005). Tru- 
dolyubov (2008) also used XMM-Newton, but considered only 
the 0.3 - 7.0 keV band pass, whilst Strohmayer & Mushotzky 
(2003) calculated the bolometric luminosity for each source. 
Finally, Liu & Mirabel (2005) collated information from pub- 
lished works and so provide details of the observed peak lu- 
minosity over an identified range specific to each source (all 
luminosities taken from this work were calculated over the 0.5 

- 10.0 keV energy range). 

We searched the Chandra archive and Hubble Legacy 
Archive (HLA) for publicly available observations of each of 
the 45 sources (using data available November 2011). Twelve 
objects which lack either Chandra or Hubble data are marked 
with a '* ' in Column 1 of Table[TJ and we do not discuss them 
further in this paper. We are left with a sample of 33 ULXs 
residing within 5 Mpc. 

3. DATA COLLATION & REDUCTION 

With the inception of the HL^] designed to optimise the 
science from HST, we are able to collate pre-processed data. 
These data sets are produced using the standard HST pipeline 
products, which combine the individual exposures using the 
IRAF task MULTIDRIZZL^\ Each field is astrometrically 
corrected (whenever possible), by matching sources in the 
field to one or more of three catalogues; Sloan Digital Sky 
Survey (SDSS), Guide Star Catalogue 2 (GSC2) and 2MASS, 
in order of preference. Information provided on the HLA 
pages states that this is only possible in ~ 80 % of ACS - 
WFC fields, due to crowding or lack of matching sources, or 
sources that are unresolved from the ground (and therefore not 
present in the catalogues). As a result, we check and improve 
on these to produce the best possible astrometry. 

To astrometrically correct these fields, it is important to 
have a number of sources in the field. We therefore opt to use 
observations with large fields of view wherever possible, and 
so select the HST instrument and detector based on these cri- 
teria. In order of preference, we select observations in avail- 
able bands from the Advanced Camera for Surveys (ACS) - 
Wide Field Camera (WFC), the Wide Field Planetary Camera 
2 (WFPC2) and ACS - High Resolution Camera (HRC). 

Another consideration is the variability that occurs within 
these systems, which would affect the emission observed in 
the optical and UV bands. From the X-ray spectra of these 
systems, we see variability on longer inter-observational time- 



http 


//heasarc.qsfc. 


nasa . gov/ Tool s/w3pimms . htmll 


http 


//hla.stsci.edu 





scales of days to years (e.g. Fabbiano 2004). Such variabil- 
ity is likely to affect the optical and UV emission, as seen in 
Galactic X-ray binary systems (e.g. Charles & Coe 2006). 
However, the X-ray variability appears to be suppressed on 
shorter, intra-observational, time-scales (Heil et al. 2009). 
Thus, we expect that near-simultaneous (< 24 hours) obser- 
vations in multiple bands are unlikely to be significantly im- 
pacted by variability. 

We fold these considerations into our observation selection 
criteria, awarding observations of multiple bands, over short 
time-scales, higher priority. Where these were not available, 
different bands were selected from multiple observations, in 
order to give a fuller view of the source. In these cases we 
must seriously consider the potential impact of optical vari- 
ability, as observations c ould have been made months, or even 
years apart (see Section |5.2| ). The observation IDs, instru- 
ment, date of observations & mode selection and filter band 
information and exposure times for selected observations are 
listed in Tabled 

We also collated and downloaded Chandra observations for 
each of our ULX sample. Wherever more than one X-ray ob- 
servation was available in the public archive, we chose the 
longest appropriate observation containing that source. In the 
case of transient ULXs this was not always the most recent ob- 
servation of the field. The observation IDs, instrument setup 
and exposure times for each of these X-ray observations are 
also listed in Table [2] 

We attempt to improve the astrometry in order to maximise 
the chance of obtaining unique optical counterparts to these 
sources. We approach the optical data first, by checking and 
improving the astrometric corrections of each field using the 
the IRAF tools CCFIND, CCMAP and CCSETWCS, in com- 
bination with either the 2MASS or USNO 2.0 catalogues (de- 
pending on number of sources available in the field). This pro- 
cess also provides the average astrometric error (3 a) across 
the field. 

We find that 14 of the 98 HST observations used in our 
analysis do not contain enough catalogued objects in the opti- 
cal/ultraviolet fields of view to allow for accurate astrometry 
corrections. In these cases we compare the field to an alterna- 
tive corrected observation in a similar waveband. We match 
sources in these observations and perform relative astrome- 
try corrections using the IRAF tools IMEXAM and GEOMAP. 
The tool GEOXYTRAN is then used to translate the position 
of the ULX to the relative field coordinates. 

In the case of NGC 4190, we find that we are unable to 
correct the astrometry of any field by matching to known cat- 
alogues. We therefore opt to take advantage of the increased 
field of view afforded us by SDSS, collecting an image of this 
region from their archive. The astrometry of this image is cor- 
rected using the 2MASS catalogue, and relative astrometry 
performed on each of the HST images. We note that where 
relative astrometry is required, the additional errors arising 
from this are also folded into our calculations. 

The astrometry of the X-ray observations must also be 
checked. We chose to use the reduced primary data provided 
by the Chandra X-Ray Center (CXC). Each observation was 
checked for any known aspect offser^J There is a small intrin- 
sic astrometric uncertainty in Chandra observations, an error 
of 0.6 "for ACIS-S, 0.8 "in ACIS-I, 0.6 "for HRC-S and 0.5 
"for HRC-I fields (90 % confidence region for absolute posi- 
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tional accuracy 1 2 1). This known error is folded into the initial 
positional error calculation. 

The X-ray astrometry can be further improved by cross- 
matching sources to the same catalogues used in the optical. 
The tool WAVDETECT, was used to identify sources in the 0.3 
-7.0 keV energy band, within 6' of the target (in some cases 
we were forced to use a smaller region, details can be found 
in Table [3]). We cross-correlated the positions of sources with 
>20 counts with the 2MASS or USNO catalogues (choice 
depending on which was used for the respective HST fields). 
Care was taken in the selection of sources, for example galac- 
tic centres were considered unsuitable as it can be difficult to 
get accurate centroiding for such sources in the optical (e.g. 
NGC 4736, Eracleous et al. 2002). 

Another concern was the limited number of sources avail- 
able for cross-matching. In cases where no suitable sources 
were found, we revert to the 90 % confidence region for abso- 
lute positional accuracy of Chandra. In some cases only one 
object was available for cross-matching in the region of the 
sky surrounding the ULX. Some previous works have used 
this single source to perform relative astrometry, but we sug- 
gest caution in doing so, as corrections can only be made 
in the XY plane, with no consideration for rotational error. 
In these cases we show both the corrected position and the 
unimproved 90 % confidence error region. All sources falling 
within the larger 90 % error region are considered, but the 
corrected position can be used to help identify the more likely 
ULX counterpart candidate. 

If multiple sources are present in both the Chandra and HST 
fields, we used a weighted average to cross-correlate the po- 
sitions of these sources, and find the required shift. When 
applying corrections with few sources, it is only possible to 
correct by shifting in the xy plane, which does not account 
for all rotational error in the telescope, but the impact is min- 
imised. As a result, errors may be underestimated. We still 
attempt such corrections but apply them with care, using them 
only as guidance in selection of 'most likely' when multiple 
counterparts lie within the larger error circle. 

When considering transient sources, the object is not al- 
ways visible in the deepest X-ray observations. In these cases 
we were able to match the position of this source to the posi- 
tion of other X-ray sources in another observation. 

Finally the accurate source positions were found using 
WAVDETECT, and any calculated shifts applied. The posi- 
tions of each ULX are listed in Table |3j along with their asso- 
ciated errors. These errors are found by combining in quadra- 
ture the astrometry errors from both fields with the source's 
individual positional error to provide the resulting error re- 
gions for each individual ULX. 

4. IDENTIFICATION OF POSSIBLE COUNTERPARTS 

The derived positions and their respective error regions are 
applied to each field in order to search for potential coun- 
terparts to these sources. Figure [T] contains a 25 x 25 arc- 
sec colour image (tri-colour wherever possible), and a finding 
chart (6x6 arcsecs) for each ULX. Error regions are plotted 
in each case, with blue representing the standard 90 % confi- 
dence ellipse, while the error derived from relative astrometry 
corrections is plotted in magenta. 

To construct the colour images, we select available wave- 
bands for each part of the optical/UV spectrum. We use fil- 
ter bands ranging from F656N to F814W to represent the red 

12 http : //cxc . harvard. edu/cal/ ASPECT /celmon/ index . h 



end of our range, filter bands F475W to F606W for the green 
band, whilst F220W - F450W are blue. In each case where 
more than one band is available, we opt for a band that gives 
the clearest view of any potential counterparts. If this is not 
required, we opt for the band with the smallest error region. 
In some cases we have no data in one or more of the colour 
bands, in which case the images are presented in one or two 
colour bands only. 

The positions of all potential counterparts are also marked 
in each field. One of these is NGC 598 ULX1, which was 
identified as the nucleus of the galaxy by Dubus et al. (2004), 
as is clearly seen in Figure p] (pg. 1). The revised Chandra 
position of IC 342 ULX2 shows that this source is also nu- 
clear. As a ULX is non-nuclear by definition, we remove both 
of these sources from our catalogue. Nine of the remaining 
ULXs have no optical counterparts in their error regions. Of 
the 22 remaining sources, 13 have a unique candidate coun- 
terpart, and the rest have up to five objects within their error 
ellipse. 

In the nine error circles which lack optical counterparts, 
limits are obtained for the approximate V band observation of 
each ULX field (listed in Table [4), given the observed back- 
ground within the positional error circle for each ULX. We 
compare these values to the expected V magnitude of each 
stellar type (Zombeck 1990) at the distance of the galaxy 
(converted to 771555 Vega-mag using SYNPHOT), with Galactic 
extinction corrections applied using the E( B — V) values in 
Table [T] O stars appear to be ruled out for 5 of these ULX. 
However, it is possible for the companion star to be a main 
sequence B star in all cases (some only valid for later-type B 
stars). We are unable to obtain a ~ V band image for M83 
XMM2, as the error region for this source lies on the ACS- 
WFC chip gap in the F555W band.We thus compare limits 
for the (only) available image, a WFPC2 F336W image, to the 
expected 771336 values for 05 V, B0 V & B5V stars (derived 
using stellar templates for the Bruzual-Person-Gunn-Stryker 
catalogue using the SYNPHOT tool CALCPHOT). O stars are 
ruled out by this comparison, but B stars are acceptable. We 
should note, however, that this does not take into account any 
extinction from the host galaxy or that is intrinsic to the sys- 
tem itself. If this extinction is high, as is seen for those sources 
in NGC 3034, this may be masking a brighter blue object. 

In those fields where potential counterparts are identified, 
we collate Vega magnitude zero points (Zpt) to allow for the 
derivation of HST filter dependent Vega magnitudes of each 
source. For observations made using the ACS instruments, 
we are able to take these values directly from Sirianni et al. 
(2005). For WFPC2 data, the HLA pipeline converts the units 
contained within the science field to electrons s _1 (like ACS) 
rather than DN (Data Number), and hence the tabulated zero 
points given in the HST Data Handbook for WFPC2 p] m ust 
also be converted, by applying a correction for the gair| i4 [ us- 
ing Zpt = tabulated zero point +2.5 x log (gain). 

Aperture photometry is performed on all potential counter- 
parts using GMApl Aperture corrections are applied to all 
fields, irrespective of instrument and detector, following the 
procedures laid down by Sirianni et al. (2005), with values 
of corrections for WFPC2 observations taken from the HST 

13 |http: //www, stsci . edu/hst/wfpc2/Wfpc2_dhb/wfpc2_ch| 
52 .html#1933986 

14 http : //www. stsci . edu/hst /wf pc2 /Wf pc2_hand_current / 
ch4 ccd!4 .html#440723 
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Figure 1. HST colour images (left) and finding charts (right) of the ULX locations. Left panel: Colour images are 25 x 25 arcseconds in size, 
over-plotted with positional error ellipse and constructed using the following filter band pass where available; blue = F220W - F450W, green 
= F475W - F606W and red = F656N - F814W. Right panel: Finding charts are made from only one individual band, showing a regions 6x6 
arcseconds, over-plotted with the combined positional error ellipse for that band (specific wave bands given in brackets after source name, 
below). Potential counterparts are highlighted numerically with associated magnitudes given in Table [4] The remainder of Figure [T]is provided 
at the end of the paper. Specific notes: displayed ULX regions are, from top to bottom, NGC 598 ULX1 (F814W) & NGC 55 ULX1 (F814W). 
NGC 598 ULX1 is contained within the nucleus of the galaxy. NGC 55 ULX1 has three potential counterparts within the error circle, however 
counterpart 3 is ruled out as it lies outside the error circle in band F606W. It will not be considered further in this paper. 

by Sirianni et al. (2005). Again, since the extinction ratios are 
also dependent on stellar type (and no blackbody is available), 
we choose to use the corrections for an 05 V star, following 
the example of Roberts et al. (2008). Although this choice 
will affect the calculated magnitudes of our sources, the im- 
pact will be minimal in most bands, with a larger (although 
still marginal) impact in the bluest bands (F435W and bluer). 
The aperture and Galactic extinction-corrected Vega magni- 
tudes are given in Table [4] 

If we wish to get a clear view of the binary system, we must 
obtain intrinsic magnitudes for these sources. To do this we 
must also take account of any absorption from either the host 
galaxy or that is intrinsic to the ULX itself. One method that 
has been used previously to correct for this extinction is to 
use the measured absorption from X-ray spectral fitting (e.g. 
Roberts et al. 2008), which will give an upper limit to the 
extinction of the optical light from the binary. 
To calculate the maximal optical extinction column, we use 



WFPC2 cookbook 

Galactic extinction corrections are also applied, using 
E(B — V) values listed in Table [T] These are used in com- 
bination with the filter specific extinction ratios depending 
on the instrument. Extinction ratios for WFPC2 data were 
taken from Schlegel, Finkbeiner & Davis (1998) where avail- 
able, or calculated using SYNPHOT when this was not possi- 
ble. The calculated extinction ratios are found by folding a 
template spectrum through the instrumental response allow- 
ing for foreground extinction using Cardelli laws (chosen for 
consistency with Schlegel et al. 1998). Although this correc- 
tion is spectrum-dependent, we find that this dependence is 
small, and we choose a 10,000 K blackbody as a first order 
estimate for these corrections, since the observed candidate 
counterparts are of unknown type. For magnitudes calculated 
using ACS fields, filter-dependent extinction ratios are given 
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the relation for X-ray-to-optical dust-to-gas ratios published 
byGuver&6zel(2009;7V# = (2.21 ±0.18) x 10 20 A v , with 
2cr errors). However, X-ray spectral fitting can be degenerate, 
and the Nh columns derived can vary substantially depend- 
ing on the author's model choice. We minimise the impact 
of model choices by using the highest quality X-ray spectra 
available, and fitting these spectra with current physical mod- 
els used to describe ULXs. We begin by collating Nh val- 
ues from published results along with their respective errors, 
preferring long observations, physical models and statistically 
good fits. When physically motivated models have not been 
applied to an object, we use values from publications of the 
deepest X-ray observation of these sources which show sta- 
tistically good fits based on phenomenological models. Our 
adopted Nh values can be found in Table [Tl Using the stan- 
dard Galactic extinction curve (Ry = AyE(B — V) = 3.1; 
Cardelli et al. 1989), we estimate the intrinsic optical mag- 
nitudes of these sources. The relevant errors are calculated 
by combining the error on Nh, obtained from literature, with 
that of Giiver & Ozel (2009). The intrinsic Vega magnitudes 
for each potential counterpart are listed in Table [4] 

We also compile details of previous identifications of poten- 
tial counterparts, in order to compare our results and search 
for the most likely candidates. Hereafter in the text, brack- 
eted values refer to the candidate counterpart ID (e.g. IC 342 
X-l (1) for candidate counterpart 1). Details of each potential 
counterpart are listed in Table|4| along with any previous iden- 
tifications of potential coun terpa rts. We compare our findings 
to previous work in Section [53 



5. COUNTERPARTS 

In our sample, we find 40 potential counterparts to the 22 
ULXs. Thirteen of these have previously been reported in the 
literature (ignoring NGC 598 ULX1), with the remaining 27 
potential counterparts identified here for the first time. Up 
to 22 of these potential counterparts may be the true coun- 
terparts, but it is possible that all the potential counterparts 
in some error circles are chance coincidences. Therefore, we 
calculate the likelihood of chance coincidences given the den- 
sity of the local stellar population of each ULX. We search for 
stellar objects within an annulus around each object with in- 
ner radius of 1" and outer radius of 3". For Circinus ULX1, 
we use a rectangle of size 3" x 8" instead, due to the chip ge- 
ometry. We expect an average of 27 ±5 objects to be present 
within the positional error regions of our ULX sample, yet 
we have observed a total of 40 potential counterparts. This 
indicates an overpopulation of 13 ±5, which is our best esti- 
mate of the number of true counterparts identified in all ULX 
fields. There is a greater likelihood of foreground object con- 
tamination when considering those ULXs in Circinus, as it 
lies in the Galactic plane, so we temporarily remove these 
sources from our calculations. If we also removed NGC 3034 
(M82) from our likelihood calculations, because the sources 
lie in an extremely obscured region of the galaxy, such that it 
is very unlikely that any optical emission would be detected 
from any real counterparts, this results in an excess of 15 ± 
4 for 36 ULXs, with 19 of these having detected candidate 
counterparts. Thus we believe that for ~ 83 % of ULXs with 
candidate counterparts, the true counterpart lies among our 
candidate counterparts. However, we note the caveat that if 
the ULX is located in a star-forming region with higher stellar 
density than areas ~ l"away, our number of true counterparts 
could be overestimated. 

In this section, we compare each candidate to stellar models 




Figure 2. Photometric stellar templates for HST filter bands, constructed us- 
ing the iraf package synphot. The templates normalised to a V band Vega 
magnitude of zero, and are grouped according to type as follows; magenta - 
O stars, blue - B-type, cyan - A type stars, green - F stars, yellow - G type, 
orange - K stars and red - M. 

to constrain the nat ure o f the donor star in these ULX binary 
In Section 



systems. 



5.1 



we consider the spectral energy dis- 
tribution (SED) and the apparent magnitude of each candidate 
to charact erise the star (assuming the light is stellar in origin). 
In Section [5T2| we check whether variability may impact our 



characterisation, which is relevant when we are in corporating 
optical data from multiple epochs. In Sectio nj5.3| we compare 
our findings to previous works. In Section [5^[ we introduce 
an additional component in the form of an accretion disc, and 
also consider the effects of irradiative heating of the star and 
disc. 

5.1. Spectral typing from stellar templates and magnitudes 

We initially consider the case where the donor star con- 
tributes 100 % of the optical light, and it has a luminosity 
and spectrum which is consistent with a single star. For these 
preliminary fits, we use only Galactic corrected magnitudes, 
ignoring any intrinsic extinction. We attempt to classify the 
donor stars by comparing the candidates with template SEDs. 
Previously, authors have done this by converting the filter 
band magnitudes to UB VRI magnitudes to compare with typ- 
ical values for different stellar types (e.g. Soria et al. 2005; 
Ramsey et al. 2006; Roberts et al. 2008; Tao et al 2011). The 
HST filters are not an exact match to other photometric sys- 
tems, which can lead to large errors in typing stars using the 
HST filter bands (see Sirianni et al. 2005 for more detailed 
discussion). 

Here, we perform typing by folding standard stellar spec- 
tra through the SYNPHOT tool CALCPHOT, a package that 
allows the user to calculate the photometric magnitudes ob- 
served for a given stellar type. In order to simplify our com- 
parisons, we choose to normalise all spectra to a V-band mag- 
nitude of zero. We use the Bruzual-Persson-Gunn-Stryker 
(BPGS) standard stellar templates associated with the syn- 
phot package. Although the atlas has the broad band cov- 
erage required for our analysis, it contains few giants/bright- 
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giants/ supergiants, which could affect our results. 

The conversion is performed for all instru- 
ment/detector/filter combinations used in the analysis of 
our ULX fields, with the resulting values plotted in Figure [2] 
Our templates range from O to M, varying in size from mam 
sequence to supergiant. The templates are grouped by colour, 
according to type (listed in figure caption). 

Simple x 2 minimisation is performed to determine the best 
fitting stellar type, whenever more than one filter band is avail- 
able, with some examples of resulting fits shown in Figure 
[3] This fitting requires that an offset be calculated in each 
instance (the shift required between its current magnitude in 
F555W and zero, for the best fitting model). This offset can 
also be considered to be the 772555 of the source (which will 
be similar in value to my). 

SED fitting is only possible when more than one filter band 
is available, however, of the 40 candidate counterparts iden- 
tified, 15 are observed in only a single filter band (due to 
chip gaps, or depth of exposures in other filters; see Ptak et 
al. 2006). Of the 25 remaining potential counterparts, we 
note that in ten instances fitting is performed where only two 
bands are available; five of our sample have three data points 
available for fitting, seven contain four HST bands and three 
contain 5 data bins for comparison to standard stellar types. 
The resultant types and offsets are displayed in Table [5j with 
sample fits shown in Figure [3] 

The best fits achieved by tnis process, given in column 3 of 
Table |5j suggest that the majority of these 25 objects are not 
best fit by OB stars, as was previously suggested (e.g. Liu et 
al. 2007, Copperwheat et al. 2007). Only two are best fit with 
O-type stars, while eight prefer B-type stars, and six are best 
fit by M stars. The best fit luminosity classes are typically not 
supergiants (only two), but main-sequence (14) or giants (9). 

We list all types that can explain the observed SEDs within 
their errors in Column 9 of Table \5\ These allow for more 
blue companions, with 20 now consistent with OB stars. Five 
of our sample appear to require later-type sources (NGC 55 
ULX1 (2), NGC 253 ULX2 (1), NGC 253 XMM6 (2), IC 342 
X-l (1) & NGC 5204 X-l (2)). However, the redder nature 
of these candidates could also be a function of reddening, as 
a result of the host galaxy or their local environment. For ex- 
ample, IC 342 X-l & NGC 3034 ULX5 have been shown to 
be heavily absorbed in X-rays on many occasions, which is 
evident from the high values of Nh listed in Table [T] Holm- 
berg II X-l may not show as high an absorption column, but 
this source resides in an excited He II region (Pakull & Miri- 
oni 2002; Kaaret et a. 2004). Although X-rays appear not 
to be heavily obscured by this nebula, it may be affecting the 
optical emission from the star. Nebulae have also been asso- 
ciated with other ULXs, including NGC 5204 X-l (Roberts et 
al. 2001) and IC 342 X-l (Pakull & Mirioni 2002). Thus, it's 
possible that the typing of these objects is incorrect, but this 
evidence for possible later-type companions to some ULXs is 
intriguing, and worth further study. 

We also consider the absolute magnitude that would be ob- 
served from each of these stellar types (Zombeck 1990; Weg- 
ner 2006). These are combined with the distance modulus 
for each source (in Table [T]) to derive the apparent magnitude 
for a star of that class at the required distanc^] We use the 

17 Exact absolute magnitudes for all star types are not given in Zombeck 
(1990) and Wegner (2006), so to compensate for this we use absolute magni- 
tudes for the most similar star type available within the text. This will lead to 
discrepancies in magnitudes of at most 1-2, depending on type. 



absolute magnitude for the V band to allow for easy compar- 
ison to our fits, and fold this through CALCPHOT to derive 
the value for the HST band F555W (i.e. 777,555 for that stellar 
type). As all templates are normalised to a V band magnitude 
of zero, this means that the choice of instrument & detector 
will have minimal impact on the derived apparent magnitude 
in this HST filter band. Since more than half of our obser- 
vations were taken by ACS using the WFC, we choose this 
combination to derive the apparent magnitudes that would be 
observed. The absolute magnitude, distance modulus and de- 
rived apparent F555W magnitudes are given in Table [5] 

As we do not have a F555W observation for every source, 
we use the offset value to derive the observed absolute F555W 
band magnitude. To calculate the value of M 555 for each can- 
didate counterparts, we combine the offset from fitting with 
the distance modulus for each potential counterpart, listing 
the resultant value in column 10 of Table \5\ These calcu- 
lations provide us with absolute magnitudes over the range 
— 1.4 < M555 < —8.2. We then used this information in 
combination with the absolute magnitudes listed in Zombeck 
(1990) to list all possible stellar types that can be observed 
at approximately this absolute magnitude (see Column 1 1 of 
Table [5]). As Zombeck (1990) only lists those classified with 
subtypes or 5 for each stellar type (e.g. 05, B0, B5, etc.), we 
classify some stars as being early (< 5) or late (> 5) within a 
stellar type. 

We compare the observed apparent magnitudes of candidate 
counterparts with the calculated magnitudes (using the offsets 
gained from x 2 fitting). To do this we group our sources into 
four categories for ease of comparison; those that differ by 
> 10 magnitudes (extreme difference), 5 < A777555 < 10 
magnitudes (large difference), those that differ by < 2 mag- 
nitudes (comparable), and those that lie in the range of 2 < 
A 777555 < 5 magnitudes (other). The ranges are designed 
to allow for magnitude and stellar template fitting errors, and 
for slight variations on tabulated absolute magnitudes while 
clearly identifying those that show striking differences. We 
find that 5 exhibit 'extreme' differences, 4 display 'large' dif- 
ferences, 9 show 'comparable' values, and 7 are classified as 
'other'. 

The sources with an 'extreme' value of A777555 are gener- 
ally best fit by later-type main sequence stars. However, in 
most cases we would be unable to see such a object, as it 
would be too faint to be observed at that distance. The fi rst is 
NGC 5204 X-l (2) and is discussed in depth in Section [5^2] 
so we will not consider it further here. The SED of NGC 253 
XMM6 (1) allows for any spectral type, while M 555 limits 
the range t o OB or later II. We will also return to this source 
in Section |5.3| The SEDs of the three remaining candidate 
counterparts only allow later-types within their errors, while 
in each case the M 555 value indicates a need for bright-giants 
or supergiants. This potential mis-classification is probably 
due to a lack of bright-giants/supergiants in the chosen cata- 
logue. 

Those sources that have a 'large' value of A777555 also ap- 
pear to be late-type sources. These sources are IC 342 X-l (1) 
& (2), Ho II X-l (1) and NGC 3034 ULX5 (1). If we again 
consider the types allowed within errors, we find that two ap- 
pear fairly well constrained (IC 342 X-l (1) and Ho II X-l 
(1)). The difference in observed and derived 777555 m ^y be 
largely due mistyping due to the catalogue used. It may also 
be the case that this mistyping is due to reddening, which must 
be considered as we are using Galactic corrected magnitudes, 
but given the increase in errors, no fits would be achievable. 
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Figure 3. Vega magnitudes of six of our sample, plotted against the stellar templates taken from the Bruzual-Person-Gun-Stryker atlas. From here it is evident 
that a range of types and data quality is available within our sample. Some are well constrained, such as Holmberg IX X-l (1), whilst others, such as NGC 3034 
ULX5 (1), have such large errors as to cover the whole range of possible stellar templates. 



Of the remaining sixteen sources, nine display similar ap- 
parent magnitudes to their offsets (Am^5 < 2) and so are 
'comparable', whilst the seven remaining sources display a 
greater divergence. In four of these cases we were unable to 
collate the exact value of My for the specified source type, 
which could induce 1-2 magnitudes of errors, so these can- 
didate counterparts can also be considered as 'comparable' 
within the increased errors. Thus, we have consistent stellar 
typing and magnitudes for thirteen of our sources. Of these 
13 sources, we find that eight are likely OB-type stars. The 
other five sources are a mixture of mid to late-type stars. The 
two potential counterparts classified as A-like stars are less 
well constrained, although in both cases any errors veer more 
towards the blue. Of the three later-type classifications (G to 
M), NGC 253 XMM6 (2) seems well constrained. This is 
very interesting, as it is unlike many previous classifications 
of ULX counterparts and as such, deserves further study. Such 
a find is also supported by the recent discovery of two near-by 
low-mass X-ray binary (LMXB) ULXs (e.g. Middleton et al. 
2012; Soria et al. 2012). Two of our ULXs (NGC 253 ULX1 
and XMM6) had multiple XMM-Newton observations in early 
2006, 8 months before the HST observations we used, which 
did not show X-ray activity. If X-ray activity did not restart 
during these 8 months, then these HST observations may be 
given an unprecedented view of the ULX companion stars, 
potentially lending more support to the LMXB hypothesis as 
other LMXB ULXs are also transients (e.g. Middleton et al. 
2012). 

Previous studies have shown that many ULXs appear to 



have blue counterparts (e.g. Roberts et al. 2008). Many of 
our potential counterparts appear to agree with this, when us- 
ing simple stellar templates. If we now consider this on a 
ULX by ULX basis (instead of each candidate counterpart 
in turn), we find that OB companions are possible in all but 
one instance - NGC 253 ULX2. This make this a key source 
for further study. While B-type companions are viable for 
all but one other ULX (NGC 5204 X-l), while we can rule 
out O-type companions in 20 cases. However, the apparent 
blue colour of these ULX optical counterparts may be due to 
the presence of a strong, blue, accretion disc component (e.g. 
Copperwheat et al. 2005; 2007). Since we did not include 
an accretion disc in these stellar template fits, finding possi- 
ble red counterparts is of interest. If these classifications are 
correct, it would suggest that the star is dominant in these fil- 
ters, as any disc emission would be intrinsically blue. This 
implies that either the ULX was dim in X-rays during these 
optical observations, or that these red objects are not the true 
counterparts to these ULXs. 

5.2. Possible optical variability; is this impacting our study? 

ULXs are observed to vary in X-rays on time-scales of days 
or more (e.g. Roberts et al. 2006), with many showing sup- 
pressed variability on time-scales of hours (e.g. Heil et al. 
2009). In Galactic X-ray binary systems, the optical emis- 
sion has also been seen to vary in a way that is related to its 
X-ray emission (e.g. Charles & Coe 2006). Investigations of 
the optical variability of these sources can be very beneficial, 
made evident by the recent work of Tao et al. (201 1), and we 



The optical counterparts of the nearest Ultraluminous X-ray sources 



9 



consider this type of analysis a valid next step for this sample, 
but it is beyond the scope of the current work. This sample 
contains only one exposure of each band for each ULX. Our 
concern here is to see if any variability in the emission from 
these sources is negatively impacting our analysis. We as- 
sume that observations taken within 24 hours should be mini- 
mally impacted, so we investigate sources that were observed 
in more than one epoch with HST. By fitting each epoch sepa- 
rately, we can investigate the potential counterpart variability 
on these time-scales, by checking for changes in their SEDs. 

Of the 31 ULX fields considered (NGC 598 ULX1 and IC 
342 ULX2 having been removed), multiple epochs were ob- 
served in 17 cases. Ignoring cases where no counterpart was 
visible, or where we only detect the candidate counterpart in 
one band, there are only eight potential counterparts to con- 
sider; NGC 4190 X-l (1), NGC 253 ULX3 (1), M81 X-6 (1) 
& NGC 1313 X-l (1), along with IC 342 X-l (1) and (2) 
and NGC 5204 X-l (1) and (2). To look for variability, we 
fit the SEDs from each epoch of observations separately with 
the same \ 2 test outlined above, noting any changes in the 
preferred stellar type. 

NGC 1313 X-l (1) is observed in five different energy 
bands, four during the same epoch. This only allows for ad- 
ditional fitting of one epoch. With the removal of the F606W 
band we find no change, suggesting little to no variation in 
the SED. A similar result is found when fitting NGC 4190 X- 
1(1). The full SED of this source contains four filters, three 
of which were scheduled together. By fitting only these three 
bands, the best fitting stellar type appears slightly redder (B9 
V). Such changes are not significant (within errors), but vari- 
ations like this have been noted before in NGC 1313 X-2 (1) 
by Mucciarelli et al. (2007) and NGC 5055 X-2 (Roberts et 
al. 2008), with variations attributed to non-stellar processes 
(possibly from the accretion disc). 

Although IC 342 X-l and NGC 253 ULX3 have multiple 
bands for each observational epoch, each of their potential 
counterparts are visible in only one band from a different 
epoch. If we remove that one band and refit, we find a sta- 
tistically similar fit, with the same range of stellar templates 
allowable (within errors). 

NGC 5204 X-l (1) and (2) have two bluer bands from one 
observation (09370_01), whilst two redder bands are from 
another (08601 _39). For each candidate counterpart we fit 
each observational epoch separately. (1) does not show any 
significant changes in fitting, while (2) shows an extreme 
change. This candidate counterpart has previously been iden- 
tified and discussed by Goad et al. (2002) as a star cluster 
(09370_01). When Liu et al. (2004; 08601_39) revisited this 
source, they incorporated higher resolution data from the ACS 
HRC, which was able to resolve the source into two compo- 
nents, revealing the presence of an 05 V star and a redder 
star cluster. We find that the complete data set is best fit by 
a template for an F8 V star, but that the complete data is not 
well described by any stellar type, with fits of the blue and red 
bands showing a two-component fit. In the blue bands we are 
seeing emission from primarily the young O-type star, whilst 
the redder bands contain emission from both the star and the 
nearby cluster. This appears to confirm the suggestions of Liu 
et al. (2004). 

Finally we return to M81 X-6. The data for its candidate 
counterpart is unlike any of our available templates. This 
spectrum appears to be bright at both the red and the blue 
ends. The initial fit for all data is a B2 III, whilst two of the 
separate observations (09073_01 & 10584.18) are fit by an A3 



III (O to F) and M0 V (O to early M) respectively. One way to 
explore this further is to split the spectrum by wave band, fit- 
ting the F336W, F435W and F555W photometric magnitudes 
in the first instance and F555W, F606W and F814W in the 
second. By doing this we note a large difference in the spec- 
tral types observed, and obtain the best fits for an 08 f and M0 
V template star respectively. This is too large a discrepancy to 
be explained by variability, and would seem to indicate some 
form of source confusion/contamination (c.f. NGC 5204 X-l 
(2)) by the combination of emission from two stars, or a two 
component spectrum that could be explained by an irradiated 
disc and a red supergiant, the second of which is an intriguing 
option. Further analysis is required to confirm either scenario. 

Of the 8 potential counterparts discussed above, we find 
minimal impact from variability in 6 cases. In the remaining 
2 cases, we find that the most extreme variations can be more 
easily explained by the presence of a two component spec- 
trum. This could encompass a star and an accretion disc, or 
it could be a the presence of multiple stars (or a star + stellar 
cluster). This shows the importance of both SED construction 
from a single observation, and for variability studies. Each 
can give us valuable information on the optical counterpart of 
the ULX, but combining multiple epochs within a single SED 
can lead to misinterpretation. 

5.3. Comparisons to previous studies 

Table [5] notes any previous identifications and source clas- 
sifications. Where more than one candidate counterpart is 
present in our sample, we list the previous identification 
alongside the counterpart matching that referred to in the lit- 
erature. Thirteen potential counterparts have been previously 
identified from our sample (discounting NGC 598 ULX1), al- 
though IC 342 X-2 (1) was unable to be classified in previous 
works. We find that we are still unable to classify it with cur- 
rent archival data. Although NGC 5204 X-l (1) & (2) and 
M81 X-6 (1) wer e pr eviously classified, and are discussed in 
detail in Section 15.21 As a result we will not discuss them 
further in this section. 

Of the nine cases remaining, we find that four of our stellar 
type ranges are in agreement with previous results (Ho IX X-l 
(IkNGC 1313 X-2 (1), NGC 2403 X-l (1) & NGC 5408 X-l 
I n two further cases, the authors attempt to take intrin- 
sic reddening into account, altering their result. This occurs 
for IC 342 X-l (1) and Holmberg II X-l (1). We find these 
sources to be of types KO IV (IC 342 X-l (1), possible types 
cover G - K range) and A3 III (Holmberg II (1), with errors 
covering B - A types), whilst previous work found these to 
be an F type supergiant (Feng & Kaaret 2008; Roberts et al. 
2008) and an OB-type possible supergiant (Kaaret et al. 2004) 
respectively. This demonstrates that intrinsic reddening can 
have a large impact in the classification of such objects. 

We are unable to classify two of the remaining cases. The 
first of these is NGC 5128 X-l, for which we detect no coun- 
terpart. The candidate counterpart to this source was initially 
identified by Ghosh et al. (2006), as an OB star. However, 
Ghosh et al. (2006) performed their own data reduction to get 
the deepest image possible. This would suggest that the data 
retrieved from the HLA is not maximised for depth of image, 
so if we wish to consider the fainter sources the data should be 

18 A recent study by Tao et al. (2011) proposes that the emission from 
some of these sources may be non-stellar in origin. Here we are only compar- 
ing stellar type classifications, we will discuss other forms of optical emission 
later in the text, hence we do not discuss their work further at this point. 
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reduced accordingly. The same problem arises with the faint 
candidate counterpart to Circinus ULX1, where we are unable 
to obtain good constraints on the magnitudes of this candidate 
counterpart. 

A recent study by Yang et al. (201 1) identified a new poten- 
tial counterpart to NGC 1313 X-l, labeled as NGC 1313 X-l 
(1) in this study. The authors studied all of the available HST 
data on this source, finding some variability in the F555W 
band on inter-observational time-scales. They attributed this 
to variations in the accretion disc, and so used only the red- 
der bands to type the companion star to this object. Their 
analysis indicated the presence of a late-type giant, possibly a 
K5-M0 II star. However, the absolute F555W magnitude ob- 
tained from this work suggests that this source is too bright 
to be explained by a star of this class. It requires a younger 
bright giant or supergiant to explain the observed luminosity 
(assuming only a stellar origin). As our initial fitting used the 
entire HST SED, our fitting is dominated by the blue compo- 
nent, which affects our classification. Their work highlights 
the need to consider the variability of these systems (which 
we tested in Section |5.2| ), and the need to consider the pres- 
ence of the accretion disc in these extreme systems. Each of 
these are discussed further in the next section of this paper. 

5.4. Accretion disc emission 

The final point that should be considered when looking at 
the magnitudes and typing of these systems is the presence of 
an accretion disc. In accreting X-ray binary systems, optical 
radiation is released from both the companion star and the ac- 
cretion disc. The presence of such a disc would increase the 
emission and change the shape of the source spectrum. The 
colour of the star and disc will also be changed by the X-ray 
irradiation of the disc and companion star. In order to explore 
this further we apply current theoretical models designed to 
describe such systems. Attempts have been made to create 
and apply such models (e.g. Copperwheat et al. 2005, 2007; 
Patruno & Zampieri 2008; Madhusudhan et al. 2008), which 
indicated that the most likely counterpart to a ULX would be 
a high-mass donor performing mass transfer via Roche lobe 
overflow, although their findings differ dramatically for the re- 
sulting black hole masses with some suggesting MsBHs while 
others prefer IMBHs. 

Copperwheat et al. (2005) used irradiation models, in com- 
bination with models of OB main sequence stars and four su- 
pergiants ranging from F to M, to explore the resulting emis- 
sion from the system. Their work indicated that the emis- 
sion from ULXs would be impacted greatly, observing a large 
brightening in the observed magnitude due to the irradiation 
of the disc and companion star (a change of ~ 0.5 - 5 magni- 
tudes, depending on the disc size, companion star type, X-ray 
hardness and the filter band). Such a change in the absolute 
magnitude of these systems could help to explain some of the 
Arasss values observed in our sample, possibly even includ- 
ing those in the Circinus galaxy. Copperwheat et al. (2007) 
applied this model to the candidate ULX optical counterparts 
known at the time, to constrain the parameters of those sys- 
tems. This assumes that we are observing a binary system 
that contains a compact object and a companion star, with the 
accretion disc being fed by Roche lobe overflow (irrespective 
of the companion star's mass). No assumption is made on the 
mass of the compact object, with masses spanning the MsBH 
and IMBH range (10 - 1000 M©), with a wide range of stel- 
lar masses and radii also available. Mass accretion rates are 
inferred from the X-ray luminosity of the system, with the op- 



tical emission incorporating light from both the irradiated star 
and the accretion disc. We refer the reader to Copperwheat et 
al. (2005; 2007) for a more detailed discussion of the model 
and its application. 

Here we apply the same model to our current sample. To 
do this, we require X-ray flux ratios for each ULX for which 
we have possible optical counterparts. Ideally this should be 
derived from data taken concurrently with that of the optical 
data. This ideal case would allow us to understand the X- 
ray emission of the system at the time our optical data were 
observed, which would have implications on the amount of X- 
ray re-ionisation. However, since we are using archival data, 
this is generally not possible, so we work on a best efforts ba- 
sis, combining the X-ray and optical data in order to obtain 
some constraints on the nature of the system. Phenomeno- 
logical models can provide general constraints on the shape 
of low quality spectra, but the absorption columns can vary 
widely depending on model choice and data quality. We use 
published results, searching for statistically sound fits to ei- 
ther Chandra or XMM-Newton data with more physically mo- 
tivated models. Whenever these are unavailable, we consider 
phenomenological fits to the data. These models are then read 
into X SPEC to derive flux ratios for each source. The model 
flux is obtained for the 0.3 - 1.0 & 1.0 - 10.0 keV ranges, 
first with the Galactic absorption column removed and then 
with the intrinsic model fit. The derived flux ratios are listed 
in Table [6] along with the relevant models and references. 

The resulting Galactic absorption corrected flux ratios are 
combined with the Galactic extinction corrected optical mag- 
nitudes to provide a multi-wavelength view of the emission 
from these systems. Likewise, the intrinsic X-ray flux ratios 
were also combined with the intrinsic optical magnitudes for 
fitting. Both Galactic corrected and intrinsic values are fit 
with models from Copperwheat et al. (2005; 2007) to con- 
sider each system in two states: one in which the donor is in 
superior conjunction and the inclination is cos(/) = 0.0 (the 
optical light is dominated by the irradiated star, there is no 
optical disc emission); and the same conjunction with an in- 
clination of cos(/) = 0.5 (both a star and disc contribution, al- 
though the ratio of these components will vary in general the 
disc contribution tends to become dominant). This is carried 
out for each potential counterpart for which we have available 
X-ray spectra and multiple optical bands. 

We have obtained magnitudes in multiple optical bands for 
25 candidate counterparts of 18 ULXs. We list constraints on 
the binary parameters for the Galactic-corrected optical mag- 
nitudes (Table |7]) and intrinsic magnitudes (Table [8}. In each 
case we consider the candidate to be the true counterpart to 
the ULX, we treated cos(/) = 0.0 and cos(/) = 0.5 cases sep- 
arately. We select some of the more well known ULXs and 
some interesting cases from both Galactic corrected and in- 
trinsic case and provide their full fits in Figure]?] In each case, 
the figure captions contain the main findings, constructed us- 
ing both the Figures and Tables [7] & [8] 

Initially we consider only the fits from the Galactic cor- 
rected optical magnitudes. As previously stated, this can be 
considered a lower limit for the extinction of these sources. 

For 10 of the potential counterparts, we are able to not only 
constrain the mass and radius of the companion at the la 
level, but also the mass of the black hole in the system, for 
certain assumed inclinations. They are NGC 55 ULX1 (1) 
& (2), NGC 253 ULX1 (3), NGC 253 ULX2 (1), NGC 253 
XMM6 (1), NGC 253 XMM6 (2), Ho IX X-l (1) & (3), IC 
342 X-l (1) and NGC 5204 X-l (2). These findings are sum- 
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Figure 4. pg. 1 - Confidence contours for select candidate ULX counterparts. Parameters are plotted in two groups of four for each ULX; 
top four containing fits achieved using Galactic corrected data while the bottom four relate to intrinsic magnitudes. Left-hand panels represent 
a superior conjunction with an inclination of cos(/) = 0.0 (observed emission dominated by irradiated star). Right hand panels are for an 
inclination of cos(7) = 0.5, again at superior conjunction (emission from both irradiated star & disc, although generally dominated by the disc). 
Panels are missing when no constraints are obtained. Here are fits for NGC 253 ULX2 (1). Galactic extinction/absorption corrected data 
suggests that Mbh > 38M© for an edge on system, where the companion mass and radius constraints indicate a late-type giant to be a likely 
counterpart (comparing to Zombeck 1990). When the inclination is set to cos(/) = 0.5, then Mbh < 59OM at the 1 a level. Here, stellar 
constraints suggest either an OB-type companion (in which case the source would need to be heavily reddened) or a later-type giant. Switching 
now to the intrinsic fits, we lose many constraints. We only obtain lower limits on the star's radius when the system is inclined, while we obtain 
lower limits on both the stellar mass and radius at cos(/) = 0.0. The mass constraints show that the system cannot be explained by a LMXB, 
but that either an intermediate or high-mass companion is possible, where stellar radius constraints tell us that we are observing either an OB 
star or a giant or supergiant (Zombeck 1990). 



marised in Table [9] When comparing these to the rest of the 
population of canaidate counterparts, we find that they gener- 
ally have smaller errors than the other potential counterparts. 
Eight of these also seem to be redder in colour, with the ex- 
ceptions being Ho IX X-l (1) & (3). 

Two of these sources, NGC 253 ULX2 (1) and NGC 5204 
X-l (2), show both upper and lower black hole mass con- 
straints, although each are provided for a different inclination. 
We find that we are only able to constrain the lower black 
hole mass limit from NGC 253 ULX2 (1) when we assume 
that all of the observed optical emission is from the irradiated 
companion star (cos(/) = 0.0). Here, fitting implies that the 
black hole must be greater than ~ 37.5 M©, a mass that is 
larger than any observed MsBH seen to date, but still cover- 



ing both the massive stellar and intermediate mass regimes. 
When we consider the alternative scenario, in which we ob- 
serve emission from both the companion star and the accre- 
tion disc (cos(z) = 0.5), we obtain an upper limit of 590 M . 
Again, this covers all classes of black hole. 

In the case of NGC 5204 X-l (2), we see the reverse. 
We obtain an upper limit on the mass when only consider- 
ing emission from the companion, while we obtain a lower 
limit when cos(/) = 0.0. These are M BH < 2OM & 
Mbh ^ 15.9M , respectively. In the first case, this would 
imply that we are observing a regular sMBH, similar to those 
seen in our own galaxy, but accreting at a much higher accre- 
tion rate (if this is the correct companion). However, if the 
inclination is increased, so that we also see some of the accre- 
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tion disc in the optical bands, the lower limit allows for any 
category of black hole. However, its two component nature is 
thought to be a result of a star cluster and an 05 star. As a 
result it is not considered the likely companion to this ULX. 

In two cases, we obtain lower-mass limits at an angle of 
cos(/) = 0.0. These are Ho IX X-l (1) and IC 342 X-l (1). 
In each case we obtain a lower limit of ~ 19 M , placing 
their lowest mass at the upper end of those observed in our 
own galaxy. In each of the remaining six cases, we obtain up- 
per black hole mass limits, under the assumption that the ob- 
served optical data incorporates both emission from the irra- 
diated disc and companion star (cos(z) = 0.5). NGC 55 ULX 
(1) & (2), NGC 253 ULX1 (3) and NGC 253 XMM6 (2) all 
have upper limits that are approximately < 400 M , such that 
no classification on black hole type can be made. However, 
NGC 253 XMM6 (2) and Ho IX X-l (3) have upper mass 
limits that lie solely within the range of massive MsBHs, with 
masses < 85M , although we should note that Ho IX X-l 
(1) is thought to be the more likely counterpart, due to He II 
emission in its optical spectra (e.g. Roberts et al. 2011). 

In cases where we have only constraints on the compan- 
ion, we find some general trends emerging. When the system 
is inclined such that the disc is edge-on (so that emission is 
purely from the companion star), we find that the mass and ra- 
dius ranges of the star tend to increase with increasing black 
hole mass. The opposite trend is present when emission is 
also thought to come from the disc. There are a few instances 
when this is not the case. In four cases (NGC 253 ULX3 (1), 
M81 X-6 (1), NGC 1313 X-l (1) & NGC 1313 X-2 (1)) we 
see the opposite trends occurring, while Ho IX X-l (1) & IC 
342 X-l (1) show approximately the same values across the 
range of black hole masses. 

What can these stellar mass ranges tell us about the system? 
Using the following approximate mass ranges, we are able to 
classify the potential companion stars of these ULXs. Low 
mass stars are considered to be those of < 1M , those in 
the range of 1 < M* < 1OM are intermediate, while those 
with M* > 1OM are considered high-mass stars. We find 
that, although it was previously thought that these systems are 
HMXB, the presence of low-mass stars cannot be ruled out in 
27 cases of the 52 tested (26 potential ULX counterparts in 2 
scenarios). Intermediate mass stars are possible in 41 cases, 
while high-mass stars cannot be ruled out in 34 cases. If we 
fold in the observed absolute magnitudes of these candidates 
(from Table [5]), we can again compare to Zombeck (1990) to 
see how many of these ULXs could be playing host to OB 
companions. Of the 18 ULXs that had potential counterparts 
available for fitting, all can hold B stars (depending on choice 
of inclination and black hole mass), while only 8 can contain 
O stars. 

Due to the increased errors on our intrinsic magnitudes, Ta- 
ble [8] shows that we are only able to obtain one black hole 
mass constraint from the 52 cases considered. The fits to Ho 
IX X-l (3) provide an upper bound of 350 M in the case of 
cos(z') = 0.5. This is also listed in Table [9] We are also un- 
able to obtain mass constraints on the companion in 30 of the 
52 cases, 14 of which also give no constraints on the radius 
of the companion. Where constraints are achieved, they fol- 
low the same trends as those outlined from Galactic extinction 
corrections, however we have many more lower limits on the 
star's mass and radius, as constraints on their upper bounds 
are lost. Where we have constraints, we attempt to classify 
the candidate counterparts as low, high or intermediate mass, 
we find that 7 cases can be described by a system containing 



a low-mass companion, that an intermediate mass star cannot 
be ruled out in 20 cases, while 22 cases may contain a high- 
mass companion. 

6. SUMMARY 

Here we present the findings of our survey of the poten- 
tial optical counterparts to ULXs, that combines data from 
both the Hubble Legacy Archive and the Chandra Space Tele- 
scope. We collate information pertaining to those ULXs re- 
siding within ~ 5 Mpc, and search for any potential coun- 
terparts. We find that from our initial sample of 45 ULXs, 
12 have no archival data. In the remaining 33 cases we col- 
lated data from each telescope and correct the astrometry of 
the downloaded data. By cross-correlating Chandra and HST 
field, we found two of the sample to reside within the nucleus 
of their host galaxies, they were therefore removed from our 
analysis. We find 22 of the 31 remaining ULXs show the pres- 
ence of candidate optical counterparts, with 13 ULXs having 
a single optical candidate in the X-ray error region. Nine of 
our sample have no observed counterpart within the error re- 
gion, although as Ptak et al. (2006) highlighted in some cases 
this will be due to insufficient depth in the exposures of these 
fields. The remaining 22 ULXs have a total of 40 potential 
counterparts, 26 of which are observed in multiple bands af- 
fording us the opportunity to attempt classification. It is ob- 
vious that not all can be the true counterparts to these ULXs, 
so we derive the number of chance coincidences to remove 
these from our sample. This suggests that 13 ± 5 of the de- 
tected counterparts are correct for the 22 ULXs considered. 
When we remove Circinus 1 and NGC 3034 sources from the 
catalogue, this changes to 15 ± 4 for 19 ULXs. 

We find that initial identifications of potential counterparts 
show no prevalence of a single stellar type. Classifications 
cover the wide range of types from blue OB stars to red M 
types, and range in size from main sequence to supergiants 
(that are possibly reddened). When considering the derived 
absolute magnitudes of these sources in the F555W filter band 
(~ My), the results are more suggestive of giants/bright- 
giants/supergiants in the majority of cases, although some ap- 
pear too bright to be explained by even the most luminous 
stars. The presence of such luminous objects indicates that, 
in some cases we are either observing foreground sources that 
are not related to the ULX, or that the stellar emission is en- 
hanced by emission from an irradiated star and/or accretion 
disc. Such emission could easily brighten the system by up to 
~ 5 magnitudes (Copperwheat et al. 2007), in agreement with 
the observed disparity in optical flux. If instead we combine 
the range of possible stellar types with the derived absolute 
magnitudes, this indicates that we are mainly observing OB- 
type stars, with OB stars ruled out for only 1 ULX - NGC 
253 ULX2. This source has only one detected counterpart, 
a red SED that can only be explained by late K or M stars, 
while its absolute V magnitude is ~ - 6.2. However, this 
magnitude and SED fitting was obtained from the Galactic- 
corrected magnitudes, so would it be reasonable in the intrin- 
sic case? If we compare galactic and total columns we find 
that N H goes from 3 to 20 (xlO 20 cm" 2 ). This means that 
E(B — V) changes from ~ 0.06 to ~ 0.8, a relatively small 
change. However, the errors on the derived E(B — V) values 
would be considerably larger due to the large uncertainties in 
Nh- This means that we are unable to explore this option at 
present. To test this further we would need great constraints 
on the extinction/absorption of this system, constraints that 
could be achieved using deeper X-ray observations. 
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The application of X-ray irradiation models provides con- 
straints on the black hole mass in only 10 cases, when fitting 
each of the potential counterparts using only the Galactic cor- 
rected X-ray and optical values (assuming in each case that 
this is the correct counterpart to the ULX). However, con- 
straints are limited. In one case, the limit suggests a sMBH 
(NGC 5204 X-l (2), cos (i) = 0.0), while in another the com- 
panion can be either a sMBH or a MsBH(Ho IX X-l (3), cos 
(i) = 0.5). These are interesting results as they agree with the 
current theory regarding these more standard ULXs. How- 
ever, in another instance, only a sMBH is ruled out (MsBH 
& IMBHs allowable; NGC 253 ULX2 (1), cos (i) = 0.0). We 
find that the fits from five of these cases provide an upper limit 
on the black hole mass of the order of hundreds of M (NGC 
55 ULX1 (1) & (2), cos (i) = 0.5; NGC 253 ULX1 (3), cos (i) 
= 0.5; NGC 253 ULX2 (1), cos (i) = 0.5) , while the three re- 
maining cases cannot rule out any classification of black hole 
(Ho IX X-l (1), cos (i) = 0.0; IC 342 X-l (1), cos (i)\ NGC 
5204 X-l (2), = 0.0, cos (i) = 0.0). We lose almost all con- 
straint in the intrinsic case, obtaining only 1 upper limit of 350 
Mq for Ho IX X-l counterpart 3 when the system is inclined. 

We also obtain companion stellar constraints in some cases 
for both galactic extinction/absorption corrected values and 
intrinsic data. We find that, although it was previously thought 
that these systems are HMXBs, the presence of low-mass stars 
cannot be ruled out in 27 cases of the 52 tested (26 potential 
ULX counterparts in 2 scenarios) for Galactic corrected val- 
ues, while 7 show that low-mass companions lie within ac- 
ceptable mass and radius ranges for the intrinsic case. Inter- 
mediate mass stars are possible in 41 cases, while high-mass 
stars cannot be ruled out in 34 cases for Galactic corrected 
magnitude/flux ratio fitting, while 20 & 33 intrinsic cases can 
be explained by intermediate or high mass stars respectively. 

This work has also highlighted several sources for which 
additional photometric or spectroscopic analysis could pro- 
vide interesting science. NGC 253 is a galaxy containing 
two transients that may have been turned off at the time of 
the archival HST observations. Another interesting thing to 
note for the companions in this galaxy, is that they appear to 
be very red, and well-fit by later- type companions. Follow- 
up photometric analysis of the stars in this galaxy could 
give greater constraints on possible companion types, while 
new deeper observations, taken with simultaneous X-ray data 
would confirm the level of X-ray emission from the transient 
ULXs, and show any change in optical emission from these 
sources. This analysis has also revealed several good can- 
didates for optical spectroscopic follow-up, 5 of which have 
been successfully awarded time with the Gemini Observatory 
as part of our ongoing programme (NGC 1313 X-2 (1), NGC 



5204 X-l (1) & Ho IX X-l (1), NGC 4395 X-l (1), & NGC 
253 ULX2 (1), a number of which will be discussed in Glad- 
stone et al. (in prep), while two others have been studied by 
alternate groups (NGC 5408 X-l (1), e.g. Cseh et al. 2011 & 
Grise et al. 2012; Ho II X-l (1) PI: Liu). 
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Figure 1. continued, pg. 2 - Specific notes: displayed ULX regions are, from top to bottom, NGC 4190 X-l (F300W), NGC 253 ULX 1 
(F475W) & NGC 253 ULX2 (F606W). 




Figure 1. continued, pg. 3 - Specific notes: displayed ULX regions are, from top to bottom, NGC 253 ULX3 (F814W), NGC 253 XMM6 
(F814W) & M81 X-6 (F435W). 




Figure 1. continued, pg. 4 - Specific notes: displayed ULX regions are, from top to bottom, Holmberg IX X-l (F435W & F330W) & NGC 4395 
ULX1 (F336W). Two finders are provided for Ho IX X-l as source 1 and 2 are only clearly separated in the F330W band due to the greater 
spacial resolution of HRC. However, source 3 is not detected in this band. 




Figure 1. continued, pg. 5 - Specific notes: displayed ULX regions are, from top to bottom, NGC 1313 X-l (F606W), NGC 1313 X-2 (F555W) 
&IC342X-1 (F555W). 
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Figure 1. continued, pg. 6 - Specific notes: displayed ULX regions are, from top to bottom, IC 342 X-2 (F435W), IC 342 ULX2 (F814W), & 
IC 342 X6 (F814W),. IC 342 X-2 has three potential counterparts within the blue error ellipse, while only two are contained within the magenta 
region. As a result we consider sources 1 & 2 to be more likely counterparts, but retain 3 in our analysis due to uncertainties on cross-matching 
fields. Cross-matching Chandra positions with HST data shows that IC 342 ULX2 is associated with the nucleus of the galaxy, and therefore 
should not be considered a ULX. It is most likely that this is a low luminosity AGN. 




Figure 1. continued, pg. 7 - Specific notes: displayed ULX regions are, from top to bottom, Circinus ULX1 (F606W), Circinus ULX3 (F606W) 
& Circinus ULX4 (F606W). 




Figure 1. continued, pg. 8 - Specific notes: displayed ULX regions are, from top to bottom, NGC 2403 X-l (F606W), NGC 5128 ULX1 
(F555W) & CXOU J132518.3-430304 (F555W). Potential counterparts 2 & 3 to NGC 2408 X-l are ruled out as they reside outside the error 
region in other bands. 




Figure 1. continued, pg. 9 - Specific notes: displayed ULX regions are, from top to bottom, NGC 4736 X-l (F814W), Holmberg II X-l (F814W) 
& M83 XMM1 (F435W). Potential counterparts 2 & 3 to M83 XMM1 are ruled out as they reside outside the error region in other bands. 
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Figure 1. continued, pg. 11 - Specific notes: 
(F814W) & NGC 3034 ULX5 (F435W). 



displayed ULX regions 



are, from top to bottom, NGC 3034 ULX3 (F555W), NGC 3034 ULX4 




Figure 1. continued, pg. 12 - Specific notes: displayed ULX regions are, from top to bottom, NGC 3034 ULX6 (F814W) & CXOU 
J095550.6+694044 (F814W). 
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Figure 4. continued, pg. 2 - Specific notes: displayed above are fits for Holmberg IX X-l (1), that was suggested by Liu et al. (2004) as the 
more likely candidate. We obtain black hole mass constraints for Galactic extinction/absorption correction data in superior conjunction with 
cos(z) = 0.0. We find that Mbh > 18.6M (assuming la constraints, see Table [7}. We find that stellar mass and radius constraints are also 
quite confined, only allowing f or an early F-type supergiant (Zombeck 1990), remarkably different from the B giant or supergiant suggested 
by our SED fitting (see Section |5Jl . The alternative is that the system is at low inclination, allowing us to see blue optical emission from the 
disc. Here, the stellar mass and radius constraints rule out the possibility of a supergiant, but still allow for a range of giants (A to mid G) and 
main sequence (B to mid G) companions. The black hole mass constraint is lost when we switch to intrinsic magnitudes, with fits ruling out 
the possibility of this system containing a main sequence star for all types except B, and the later- type supergiants (G or later), for the case of 
cos(/) = 0.0. At cos(/) = 0.5, we see that the F-type supergiants are also ruled out, but that the entire main sequence is available, depending on 
the assumed black hole mass. 
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Figure 4. continued, pg. 3 - Specific notes: displayed above are fits for Holmberg IX X-l (3). Here we see that the star's radius allows us to 
place constraints on the black hole mass in all cases, with upper limits ranging from 85 - 350 M (at 1 a level; see Table [7] & [8}. If we assume 
that the top four panels are correct (Galactic extinction correction only), we rule out IMBHs, and the possibility of an HMXB. This is because 
the stellar constraints are M* < 6.2M and i?* < 2.5Rq (for Mbh — 10M©). When we switch to the intrinsic scenario (no shielding of the 
star), the black hole mass constraints relax so that both MsBHs & IMBHs are possible, but the star's constraints are still kept to M* < 12.6M 
and R* < A.2Rq. The companion could be a mid-B or later main sequence star, or a mid-class (~ F) giant. 
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Figure 4. continued, pg. 4 - Specific notes: fits for NGC 1313 X-2 (1). No constraints are found for the black hole mass in this ULX, but 
some constraints are achieved for the companion. Also, if we incorporate findings from studies of the surrounding stellar population, we can 
use an upper mass limit of ~ 12M (Grise et al. 2008). This suggests a companion mass range of 1.7 < M* < 12M and a radius range 
of 2.6 < R* < 325i?0, when using galactic extinction / absorption corrections and cos(/) = 0.0. This allows for early B to F0 type main 
sequence stars, along with A to mid-F giants (III) (Zombeck 1990). If we switch to the alternative inclination, we find that the range of stellar 
radii is reduced while lower limits on the mass range are lost. This implies that early-type main sequence stars are permissible, along with A to 
early G-like giants, while all supergiants are ruled out. Intrinsic data suggests only a lower-mass limit for the companion in the case of cos(/) = 
0.0, while the radius range is large. Combining these with Grise et al. (2008), we rule out all main sequence stars and all but F type supergiants, 
however, any giants are acceptable (using Zombeck 1990). If we switch to the inclination of cos(/) = 0.5, we lose stellar mass constraints but 
retain constraints on the star's radius, allowing any main sequence stars, all giants, but only OB-type super giants. Again, combining this with 
the mass limit 12 M©, this rules out O and early B stars of all classes. 
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Figure 4. continued, pg. 5 - Specific notes: displayed above are fits for IC 342 X-l (1). Fitting provides a lower limit on the black holes 
mass at the 1 a level in the upper two left hand plots (galactic corrected with inclination of cos(/) = 0.0. This limit indicates a lower-mass of 
Mbh > 18.6M©. The star mass and radius ranges, provided in the top four panels, are also narrow enough that they can also place good 
constraints on the companion star type. It indicates that we are most likely observing a late G or early K supergiant when cos(/) = 0.0 and 
a sliding scale for cos(/) = 0.5, covering a wide range of stellar types. For sMBHs the donor must be a high-mass companion that is either 
giant or supergiant in class. However, an M555 = —5.7 (from Table [5}, rules out type la stars, although many lb's are still allowable. If 
Mbh = 100M Q the companion can be an O or early B main sequence, or any giant (class III) star. Finally, if we are observing an IMBH 
(Mbh = 1OOOM ), the companion can be an OB main sequence, or a mid to late B to G III. Now we turn to the intrinsic magnitude for this 
source at an inclination of cos(/) = 0.5, where we are unable to obtain mass constraints once again. We also find that the radius range generally 
increases with increasing black hole mass, allowing for O & early B main sequence stars, any III, or M0 and younger I, with the smallest range 
narrowing this to O-like main sequence, B or K III, K0 or younger supergiants. 
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Figure 4. continued, pg. 6 - Specific notes: displayed above are fits for IC 342 X-l (2). By combining these plots with the values from Tables 
[7]&|8] we obtain stellar mass constraints of 1.3 < M* < 51.3 M for Galactic corrected magnitudes with radius constraints 3.5 < R* < 
TJI.oRq for an edge on system, which drops to M* < 1.5M© and 5.1 < R* < 641.6 Rq when intrinsic magnitudes are used. This allows 
for OB main sequence, any III, or B to mid K supergiant stars in the first instance, and O or early B main sequence, any III, or all but the 
reddest supergiant stars in the second instance. If we instead switch to cos(/) = 0.5, the star's mass and radius constraints provide us with a 
high upper limit on the mass of M* < 44. 9M© with a lower radius limit of i?* < 12. 9Rq, with the radius decreasing with increasing black 
hole mass. This combination rules out all O classifications, all supergiants, and K or early B III. The intrinsic magnitudes provide us with even 
less constraints, obtaining only an upper limit on the companion star's radius (i?* < 55. 5i?©) allowing for any main sequence star, any giant, 
or mid A or younger supergiants. 
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Galactic corrected parameters 
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Figure 4. continued, pg. 7 - Specific notes: displayed above are fits for NGC 5204 X-l (1). We obtain mass and radius constraints for three 
of the four scenarios considered, losing mass constraints only in the intrinsic inclined (cos(7) = 0.5) case. For cos(/) = 0.0, in the galactic 
extinction corrected case the mass and radius constraints listed in Table [7] constrain the allowable companion star types to O V, OB III or early 
B I, while intrinsic magnitudes extend the range to O or early B main sequence stars, any giants, or early K or younger supergiants. When we 
switch to an inclined system, we once again obtain a sliding radius scale for the companion star, which decreases with increasing black hole 
mass. For lower-mass black holes, we find strong constraints on the possible companion star, only allowing for early B supergiants. As the 
black hole mass increases, however, we find that none are comparable with the mean values (Zombeck 1990). This would indicate that either 
the the companion is not of a mean classification or that this is not the correct scenario or counterpart for this ULX. This could also be due to 
the fact that the assumption of galactic corrected magnitudes is incorrect in this case. The intrinsic case, however, allows little to no constraints, 
meaning that any companion type is possible. 
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Table 1 

The ULX sample listed by distance. 



Source 


Alternative names 


RA 
(J2000) 


Dec 
(J2000) 


d a 
(Mpc) 




E(B-V) C 




(10 39 er 


e 

gs- 1 ) 


NGC 598 ULXl^) 


Source 3^ 7 ^ 
Source l^ 10 ^ 
M33-X-8< JJ ) 

CXOUJ013351.0+303937( iS ) 


01 33 50.9 


+30 39 37 


r\ ft} 


4.59 


0.042 


3.6 E( y ) 


2.5 


a:XL 7 J 


NGC 55 ULXl(^) 


NGC 55 6^3) 

XMMU J001528.9-391319(^) 


00 15 28.9 


-39 13 19 


1.74(2) 


1.71 


0.014 


23.9 E( p ) 


1.3 




NGC 4190 Xl( 5 ) 




12 13 45.4 


+36 37 55 


2.8( 5 ) 


1.54 


0.030 


12T (i5) 


2.31 


/3RP(3) 


NGC 253 ULXl^) 


NGC 253 PSX-l^) 
CXOUJ004734.0-251637( ie ) 


00 47 34.0 


-25 16 37 


3.0^ 


1.42 


0.019 


N/A 


1.2 




NGC 253 ULX2^) 


NGC 253 PSX-2^) 
NGC 253 X2( 5 ) 
NGC 253 XMMl^) 
CXOUJ004733.0 
-251749( ie ) 
S102( i ) 


00 47 33.0 


-25 17 49 


3.0( 5 ) 


1.42 


0.019 


20 TU) 


5.7 


aXlU) 


NGC 253 ULX3^) 
NGC 253 X20( 5 )* 


NGC 253 PSX-3^) 
NGC 253 XI ( 5 ) 
CXOUJ004733.4-251722( ie ) 
NGC 253 ULXl( 5 ) 


00 47 33.4 
00 48 20.0 


-25 17 22 
-25 10 10 


3.0( 5 ) 
3.0( 5 ) 


1.42 
1.52 


0.019 
0.019 


N/A 
N/A 


1.74 
2.37 


/5RP(3) 
/3RP(3) 


NGC 253 XMM2^>* 


NGC 253 X9( 5 ) 


00 47 22.4 


-25 20 55 


3.0(3) 


1.41 


0.019 


39 T(^) 


2.7 


aXlU) 


NGC 253 XMM4^>* 




00 47 23.3 


-25 19 07 


3.0(3) 


1.42 


0.019 


1.2 T^) 


2.2 


aXl(^) 


NGC 253 XMM5^>* 


NGC 253 X7( 5 ) 


00 47 17.6 


-25 18 12 


3.0(3) 


1.42 


0.019 


3.4 tU) 


2.2 


aXlU) 


NGC 253 XMM6^) 


NGC 253 X6( 5 ) 

RX J004742.5-25150l( i7 ) 


00 47 42.8 


-25 15 06 


3.0(3) 


1.43 


0.019 


3.9 T(^) 


3.1 


aXl(^) 


M81-X-6^) 


NGC 3031 X9( 5 ) 
NGC 3031 ULXl( 5 ) 
M81 XMMl^) 
Source 7 ( 7 ) 
CXOUJ095532.98 
+690033.4^ /; 


09 55 33.0 


+69 00 33 


3.4(3) 


4.16 


0.080 


19 


3.84 




Hoi IX X-l(^) 


NGC 3031 ULX 2^) 

M81-X-9^) 

HolIXXMMl^) 

Source 17( 7 ) 

NGC 3031 I0^ loj 

H44( i5 ) 

IXO 34( 20 ^ 


09 57 54.1 


+69 03 47 


3.42(3) 


4.06 


0.079 


12.1 


13.4 


7 RP(3) 


NGC 4395 ULXl( J > 


NGC 4395 X-l( 5 ) 
NGC 4395 XMMl^ 4 ' 
Source 12( 7 ) 
NGC 4395 X2( i5 ) 
IXO 53<* J > 


12 26 01.9 


+33 31 31 


3.6(3) 


1.35 


0.017 


10.0 e( 22 ) 


1.73 


7 RP(3) 


NGC 1313 X-l^) 


NGC 1313X2( 5 ) 
NGC l3l3ULXl( 5 ) 
Source 4( 7 ) 
IXO 7^) 


03 18 20.0 


-66 29 11 


3.7(3) 


3.90 


0.110 


21.1 E(^) 


1.3 


aXlU) 


NGC 1313 X-2^) 


NGC 1313ULX2^) 
NGC 1313X7( 5 ) 
NGC 1313ULX3( 5 ) 
Source 5 ( 7 ) 
IXO 8^ 


03 18 22.3 


-66 36 04 


3.7(3) 


3.90 


0.085 


21 EW 


4.2 
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Source 


Alternative names 


RA 
(J2000) 


Dec 
(J2000) 


d a 
(Mpc) 




E(B-V) C 




(10 39 ergs" 1 ) 


XMM J031747.6 
-663010( 5 )* 




03 17 47.6 


-66 30 10 


3.7(5) 


3.83 


0.109 


23T(5) 


1.6 


eX(5) 


I O/lOV 1(6} 

IC 342 X- 1A ; 


PCjC138zo-ao w; 
PGC 13826 ULX3( 5 ) 
IC 342 XMMl^) 
Source 2( 7 ) 
IXO 22^) 
CXOUJ034555.7 

, /TOA/) erf 23") 

+680455 ^ ; 
Source 19^) 


no ac cc in 

03 45 55.17 


i £.Q A/I CO 

+68 04 58.6 


a a(3) 


31.1 


0.565 


cc n(9) 
55 cA > 


12.77 


7 R P (5) 


IC 342 


PGC13826-X7( 5 ) 
IC 342 AJV1JV12^ 4; 
IC 342 

Source 25^) 
IC 342 X-13^5) 


03 46 15.0 


+68 11 11.2 


3.9< 3 > 


29.7 


0.559 


38T(^) 


8.43 


aXl(^) 


IC 342 ULX2^) 


PGC 13826 XI ( 5 ) 
IC 342 XMM3^> 
Source 38^4) 


03 46 48.6 


+68 05 51.0 


3.9< s > 


30.2 


0.558 


97 jU) 


2.57 


7 R P (5) 


IC 34z A-4 V y 


FCC 1 3 ozo- A3 v 1 
PGC 13826 ULX2( 5 ) 
Source 36^) 


AO /l/T AC. CA 

03 4o 4j.j4 


+05 09 51. / 


3.y y J 


on c 


n cct 
0.5 5 / 


A a t(24) 


1.49 


7 RP(5) 


IC 342 X-6( e ) 


PGC13826-X2( 5 ) 
PGC 13826 ULXl( 5 ) 
IC 342 XMM4^) 
Source 44^) 


03 46 57.17 


+68 06 22.4 


3.9( 5 ) 


30.0 


0.558 


21T(^) 


1.21 


7 R P (5) 


Circinus lULXl^) 


CG-X-l^) 

ULX42(^) 

CXOUJ141312.3 

-652013^ 7 ) 
U43^) 


14 13 12.3 


-65 20 13.0 




55.6 


1.488 


24 e( 25 ) 


24 3.72 


5C( 7 ) 


CircinusULXS^) 


CXOUJ141310.3 

-652017( S7 ) 


14 13 10.3 


-65 20 17.0 


4 (5) 


55.6 


1.468 


<487 T^ 7 ) 


1.4 


paw 


CircinusULX4( i ) 


CXOUJ141310.4 

-652022^ 7 ) 


14 13 10.4 


-65 20 22.0 


4(5) 


55.7 


1.464 


91 j(27) 


2.09 




CircinusXMMl^)* 


Source l( 7 ) 


14 12 54.2 


-65 22 55.3 


4 (5) 


55.2 


1.028 


101 t(^) 


23 


aXlU) 


CircinusXMM2^)* 




14 12 39.2 


-65 23 34.3 


4(5) 


55.1 


0.891 


112 T^) 


10.7 


aXl(^) 


CircinusXMM3^)* 




14 13 28.3 


-65 18 08.3 


4 (5) 


56.0 


1.363 


135 T(^) 


14.5 


aXl(^) 


NGC 2403 X-l( J ) 


NGC 2403 X2( 5 ) 
NGC 2403 ULXl( 5 ) 
NGC 2403 XMMl^) 
Source 6^ 
Source 21 ( 29 ) 
CXOUJ073625.5 
+653540( 50 ) 


07 36 25.55 


+65 35 40.0 


4.2( s ) 


4.17 


0.040 


23.4 E(^) 


1.73 




NGC 2403 XMM2^>* 


CXOUJ073650.0 

+653603( 5i ) 


07 36 50.2 


+65 36 02.1 


4.2( s ) 


4.13 


0.040 


18 T^) 


1.6 


aXl(^) 


NGC 5128 ULXl( 5 ) 


NGC 5128 X4( 5 ) 


13 25 19.9 


-43 03 17.0 


4.21( 5 ) 


8.37 


0.115 




9.27 


7 Rp(^) 



IXO-76^) 
ULX 40^) 
U35^) 

CXOUJ132519.9 
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Source 


Alternative names 


RA 
(J2000) 


Dec 
(J2000) 


d a 
(Mpc) 


Nv, b 




Ntt d 


(10 39 ergs" 1 ) 


NGC5128 ULXl( 5 ) 


-430317( 5S ) 


















NGC5128 X37( 5 )* 




13 26 26.16 


-43 17 15.6 


4.21( 5 ) 


8.71 


0.109 


N/A 


1.96 


7 RP(5) 


NGC5128 X38( 5 )* 




13 26 56.81 


-42 49 53.6 


4.2l( 5 ) 


8.41 


0.095 


N/A 


1.97 


7 RP(5) 


CXOUJ132518.3 
-430303 (*) 




13 25 18.3 


-43 03 03 


4.21(5) 


8.63 


0.110 


10E(^) 


1.4 


aC^ 


NGC 4736 XMMl^) 


NGC 4736X-l( e ) 
NGC 4736 X-4( 55 ) 
CXOUJ125050.3 
+410712(54) 


12 50 50.2 


+41 07 12.0 


4.3U) 


1.44 


0.018 


20T (^) 


17.9 


7XIU) 


Holmbergll X-l^^) 


PGC 23324 ULXl( 5 ) 
Hoi II XMMl^) 
Source 28 ( 7 ) 
IXO-3l(^) 
CXOUJ081928.99 
+704219.4(5°) 


08 19 30.2 


+70 42 18.0 


4.5( s ) 


3.41 


0.032 


7.9 E (9) 


17 




M83 XMMl^) 


NGC 5236 XI l( 5 ) 
NGC 5236ULXl( 5 ) 
Source 24 ( 7 ) 
IXO-82^) 
CXOUJ133719.8 
-295349^ ) 


13 37 19.8 


-29 53 49.8 


4.7(5) 


3.78 


0.066 


6.5 E(^) 


2.8 


aXl(^) 


M83 XMM2^> 


NGC 5236X6(5) 
CXOUJ133659.5 
-294959^) 


13 36 59.4 


-29 49 57.2 


4.7(5) 


3.96 


0.066 


N/A 


3.4 


«xi( 5 ) 


NGC 5204 X-l^) 


Source 23 ( 7 ) 
IXO 77^) 
CXOUJ132938.61 
+582505.6(5°) 


13 29 38.6 


+58 25 06.0 


4.8( s ) 


1.39 


0.013 


3.6 E( fl ) 


4.4 


aX(«) 


NGC 5408 X-l(^) 


NGC 5408 ULXl( J > 
NGC 5408 XMMl^) 
Source 25 ( 7 ) 
J140319.606 
-412259.572( 55 ) 


14 03 19.61 


-41 22 59.6 


4.8U) 


5.67 


0.069 


2.9 E( fl ) 


10.9 


aXlU) 


NGC 3034 ULX3^) 
NGC 3034 ULX4^) 


CXOUJ095551.2 
+694044(5°) 
ULX 14^ 5 ) 
CXOUJ095551.07 
+694045(5°) 


09 55 51.4 
09 55 51.1 


+69 40 44.0 
+69 40 45.0 


5.2< 3 > 
5.2( s ) 


3.98 
3.98 


0.159 
0.159 


1300 T( 5 °) 

320 j(36) 


12 
17 




NGC 3034 ULX5^) 


ULX 13^°) 
M82 X-l< 57 ) 


09 55 50.2 


+69 40 47.0 


5.2< s > 


3.98 


0.159 


U3T (38) 


100 


Cx(5^) 


NGC 3034 ULX6^) 


NGC 3034 XI ( 5 ) 


09 55 47.5 


+69 40 36.0 


5.K S ) 


3.99 


0.160 


61 A T(5) 


11.6 


/3CI(5£) 



CXOUJ095546.6 
694037(5^) 
CXOM82 J095547.5 

694036(5^) 



CXOUJ095550.6 
+694044(5^) 



0.55 50.6 +69 40 44 5.2( 5 ) 3.98 0.159 260 T( 57 ) 11.6 /3Cl(5 7 ) 
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Source 



Alternative names 



RA 
(J2000) 



Dec 
(J2000) 



d a 
(Mpc) 



N H b E(B-V) C 



(10 39 ergs" 1 ) 



Note. — "Distance to host galaxy collated from literature search, with cut applied at 5 Mpc. NGC 3034 has been included due to 
uncertainties in this measurement. b Galactic absorption column (in units of 10 20 cm -2 ) from Dickey & Lockman (1990) using the NASA 
HEASARC tool (this can be found at heasarc.gsfc.nasa.gov/cgi-bin/Tools/w3nh/w3nh.pl) and c Galactic extinction values found using Schlegel 
et al. (1998) via NED extinction calculator (located at ned.ipac.caltech.edu/forms/calculator.html). Each are calculated at the position of the 
ULX. d Total (T) or extra-galactic (E) absorption column (in units of 10 20 cm -2 ) found via literature search. In cases where no value was 
found we list it as not available (N/A). e Luminosity of ULX, energy band, telescope & comments relating to its calculation. Energy bands 
are noted as follows; a = 0.3 - 10.0 keV, = 0.5 - 10.0 keV, 7 = 0.3 - 8.0 keV, 6 = 0.5- 8.0 keV, e = 0.3 - 7.0 keV & C = bolometric 
luminosity. Telescope notation; C - Chandra, R - ROSAT, X - XMM-Newton. Luminosity annotation: I - intrinsic/unabsorbed luminosity, P 
derived using WEBPIMMS. * Removed from sample due to lack of Chandra and/or HST data, see Section|3]for further details. Numbers shown 
in superscript relate to the following references: (jZ) Liu & Mirabel (2005), (s) Stobbart et al. (2006), (5) Liu & Bregman (2005), (4) Winter et al. 
(2006), (5) Trudolyubov (2008), (e) Roberts & Warwick (2000), (7) Feng & Kaaret (2005), (5) Sivakoff et al. (2008), (9) Gladstone et al. (2009), 
(i0) Schlegeletal. (1998), (ii) Trinchieri et al. (1988), Grimm et al. (2005), (i5) Read et al. (1997), (i4) Stobbart et al. (2004), (i5) Roberts 
et al. (in prep), (^Humphrey et al. (2003), (i7) Barnard et al. (2008), (i5) Vogler & Pietsch(1999), (i9) Radecke (1997), (S0) Immler & Wang 
(2001), ^Colbert &Ptak (2002), ( ^ } Stobbart et al. (2006), ( 23 ) Roberts et al. (2004), ( ^ } Maketal. (2011), (S5) Kong (2003), ( ^ } Ptaketal. 
(2006), (S7) Bauer et al. (2001), ( ^ } Bergriea et al. (2008), ^Schlegel & Pannuti (2003), (50) Swartz et al. (2004), (5i) Yukita et al. (2007), 
(5S) Kraft et al. (2001), (55) Eracleous et al. (2002), ^Colbert et al. (2004), (55) Kaaret et al. (2003), (5e) Kong et al. (2007), (57) Strohmayer 
& Mushotzky (2003), (55) Kaaret et al. (2006), (59) Griffiths et al. (2000). 
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Table 2 

Details of observations used to locate possible optical counterparts to our ULX 

sample. 



Source 


Chandra 


Expos, (ks) 


HLA 


Date of 


Instrument 


Filter 


Expos. (s) 




Obs ID 


Obs ID 


Instr. 


HST obs 






NGC 598 TJT X1 


6376 


93 


06038 02 

vjvjvj^j o _u z. 


1995-10-02 


WFPC2 


F170W 


1800 






ACIS-I 


06038_02 


1995-10-02 


WFPC2 


F336W 


1800 








06038_02 


1995-10-02 


WFPC2 


F439W 


600 








06038 02 

vjvjvj^j o _vj z. 


1995-10-02 

y y y ^) y VJ UZ 


WFPC2 


F555W 


160 








05464 05 

VJ ^JtVJt —VJ^J 


1 994-09-26 

y y y^ vj y z^vj 


WFPC2 


F814W 


1280 


NGC 55 TJT X1 


2255 


59 


09765 03 

VJ J 1 UJ _Vy ~J 


2003-09-23 

L^VJVJ VJ y Z* *J 


AfS-WFr 


F606W 


400 






ACIS-I 


09765_03 


2003-09-23 


ACS-WFC 


F814W 


676 


NGC 4190 X-1 

i > vj v_x t" j — /vj y v ± 


8212 


25 


1 1012 02 

1 1 U 1 —VJ Z-! 


2008-01 -07 

jLsVJVJ O Vy 1 VJ / 


WFPC2 


F300W 


4400 






HRC-I 


1 1012 02 

1 1 Vy 1 Z- —VJ 


2008-01 -07 

jL^VJVJ O Vy 1 Vy / 


WFPC2 


F450W 


4400 








11012.02 


2008-01-08 


WFPC2 


F606W 


1600 








10905_04 


2008-03-21 


WFPC2 


F814W 


4400 


NGC 253 ULX I 


969 


14 


0521 1_01 


1994-05-29 

i y y^ vj^j y 


WFPC2 


F336W 


820 






ACIS-S 


10915 98 


2006-09-13 

L^VJVJVJ vj y X *J 


ACS— WFC 

rvv_/j vv x v^ 


F475W 


1482 








0521 1 01 

VJ^J Z- 1 1 _U X 


1994-05-29 

1 y/T Vy^y Z* V 


WFPC2 


F555W 


820 








10915 98 


2006-09-1 3 

ZwVyVyVJ VJ y X ^y 


AfS-WFr 

AV_/J VV X 


F606W 


1508 








10915_98 


2006-09-13 


ACS-WFC 


F814W 


1534 


NGC 253 TJT X2 

1>VJV Z*,J ~J 


3931 


82 


10915_98 


2006-09-1 3 

jLsVJVJVJ VJ y X *J 


AfS-WFr 

AV_/J VV X 


F475W 


1482 






ACIS-S 


10915 98 

i U7 i j — y o 


2006-09-13 

^VJVJVJ vj y X *J 


ACS-WFC 

rvv_/j vv x 


F606W 


1508 








10915_98 


2006-09-13 


ACS-WFC 


F814W 


1534 


NGC 253 ULX3 


3931 


82 


0521 1_01 


1994-05-29 

i y y^ vj^j y 


WFPC2 


F336W 


820 






ACIS-S 


10915 98 

1 U7 1 J — y O 


2006-09-13 

l^vjvjvj vj y x ~j 


ACS-WFC 


F475W 


1482 








06440 01 

VJVJ 1 1 VJ —VJ 1 


1 997-07-09 

X y y 1 VJ 1 VJ y 


WFPC2 


F502N 


2400 








0521 1 01 

Vy^y Z- 1 1 _VJ X 


1994-05-29 

1 y/T VJ^J y 


WFPC2 


F555W 


820 








10915_98 


2006-09-13 


ACS-WFC 


F606W 


1508 








10915_98 


2006-09-13 


ACS-WFC 


F814W 


1534 


NGC 253 XMM6 


3931 


82 


10915_97 


2006-09-09 

l^vjvjvj vj y vj y 


ACS-WFC 

iVVvJ VV X 


F475W 


1482 






ACIS-S 


10915_97 


2006-09-09 

^vjvjvj vj y vj y 


ACS-WFC 

rvv_/j vv x 


F606W 


1508 








10915_97 


2006-09-09 


ACS-WFC 


F814W 


1534 


M8I X-6 


735 


50 


06139 01 

vjvj x *j y —Vj x 


1995-01-31 

X y y *J VJ 1 ^y X 


WFPC2 


F336W 


1160 






ACIS-S 


10584 18 

X VJU o 


2006-03-22 

^LVJVJVJ VJ^y 


ACS-WFC 

AV_/J VV X V_x 


F435W 


1200 








09073 01 

VJ yVJ 1 —VJ X 


2001-06-04 

Z-Vy Vy X VJVJ VJ'-x 


WFPC2 


F555W 


8000 








10584.18 


2006-03-22 


ACS-WFC 


F606W 


1200 








09073_01 


2001-06-04 


WFPC2 


F814W 


8000 


Hoi IX X-1 


9540 


25 


09796 03 

vj y I yvj—VJ^j 


2004-02-07 

Z*VJVJt- Vy Z- Vy / 


ACS— HRC 


F330W 


2760 






ACIS-S 


09796 03 

VJ y 1 yVJ—VJ^j 


2004-02-07 

Z-Vy Vy^ Vy Z- VJ 1 


ACS-WFC 

iiVvJ VV X V_x 


F435W 


2520 








09796 03 

VJ y 1 y VJ —VJ ^7 


2004-02-07 

Z-Vy Vy^ Vy Z- Vy / 


ACS-WFC 

AV_/J VV X V_x 


F555W 


1160 








09796_03 


2004-02-07 


ACS-WFC 


F814W 


1160 


NGC 4395 TJT XI 

liVJV^ ^^Jy*J U J-//V1 


402 


1.2 


09774 ab 

VJ y 1 1 '-x —<X\J 


2004-06-12 

Z-Vy Vy^ VJVJ X Z* 


WFPC2 


F336W 


2400 






ACIS-S 


09774 0b 

VJ y 1 1 t- _VJ U 


2006-06-12 

Z^VJVJ VJ VJVJ X z* 


ACS-WFC 

AV_/J VV X V_x 


F435W 


680 








09774 0b 

vj y i i t- _vy u 


2006-06-12 

ZrfVJVJVJ VJVJ X z^ 


ACS-WFC 


F555W 


680 








09774_0b 


2006-06-12 


ACS-WFC 


F814W 


430 


NGC 1313 X-1 


2950 


20 


09796 a1 

VJ y 1 yvj —CI x 


2003-1 1-17 

ZrfVJVJ^y XX X / 


ACS— HRC 


F330W 


2760 






ACIS-S 


09774 05 

vj y i i 


2004-07-17 

Z^VJVJT^ VJ / X / 


ACS-WFC 


F435W 


680 








09774 05 

VJ y 1 l '-r_VJU 


2004-07-17 

Z*VJVJ'-T VJ 1 LI 


ACS-WFC 

AV_/J VV X V_x 


F555W 


680 








10210 06 

X l/Z X VJ _\J VJ 


2004-10-30 

ZrfVJVjT^ X VJ *JVJ 


ACS-WFC 

AV_/J VV X V_x 


F606W 


1062 








09774 05 

vj y i i t"_vj,j 


2004-07-17 

ZrfVJVjT" VJ / X / 


ACS-WFC 

iiVvJ VV X 


F814W 


676 


NGC 1313 X-2 


3550 


14 


09796_a2 


2003-11-22 


ACS -HRC 


F330W 


2760 






ACIS-I 


09796 02 

vj y I yvjsjL* 


2003-11-22 


ACS-WFC 

rivj vv x 


F435W 


2520 








09796 02 

VJ y 1 yVJ—VJZ* 


2003-11-22 


ACS-WFC 

iVVvJ VV X 


F555W 


1160 








09796_02 


2003-11-22 


ACS-WFC 


F814W 


1160 


IC 342 X-1 


7069 


58 


10579_b3 


2005-09-02 


ACS -HRC 


F330W 


2900 






ACIS-S 


10768_02 


2005-12-16 


ACS-WFC 


F435W 


1800 








10768_02 


2005-12-16 


ACS-WFC 


F555W 


1080 








10579_13 


2005-09-02 


ACS-WFC 


F606W 


1248 








10768_02 


2005-12-16 


ACS-WFC 


F658M 


2400 








10768_02 


2005-12-16 


ACS-WFC 


F814W 


1080 


IC 342 X-2 


2936 


2.8 


10579_a5 


2005-09-02 


ACS -HRC 


F330W 


2900 






HRC-I 


10579_15 


2005-09-02 


ACS-WFC 


F435W 


1248 








10579_15 


2005-09-02 


ACS-WFC 


F606W 


1248 
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Source 


Chandra 


Expos, (ks) 


HLA 


Date of 


Instrument 


Filter 


Expos. (s) 




Obs ID 


Obs ID 


Instr. 


HST obs 








IC 342 ULX2 


7069 


58 


06367_03 


1996-01-07 


WFPC2 


F555W 


520 






ACIS-S 


05446_0j 


1994-11-26 


WFPC2 


F606W 


160 








06367_03 


1996-01-07 


WFPC2 


F656N 


800 








06367_03 


1996-01-07 


WFPC2 


F675W 


60 








06367_03 


1996-01-07 


WFPC2 


F814W 


520 


IC 342 X-6 


7069 


58 


08199_01 


1999-08-14 


WFPC2 


F555W 


2600 






ACIS-S 


05446_0j 


1996-01-07 


WFPC2 


F606W 


160 








08199_01 


1999-08-13 


WFPC2 


F814W 


2600 


Circinus ULX1 


356 


25 


07273_01 


1999-04-10 


WFPC2 


F502N 


1800 






ACIS-S 


07273 _01 


1999-04-10 


WFPC2 


F547M 


60 








06359_08 


1996-08-11 


WFPC2 


F606W 


600 








07273_01 


1999-04-10 


WFPC2 


F656N 


1600 








07273_01 


1999-04-10 


WFPC2 


F814W 


40 


Circinus ULX3 


356 


25 


09379_64 


2002-12-11 


ACS-HRC 


F330W 


1200 






ACIS-S 


07273_01 


1999-04-10 


WFPC2 


F502N 


1800 








07273 _01 


1999-04-10 


WFPC2 


F547M 


60 








06359_08 


1996-08-11 


WFPC2 


F606W 


600 








07273_01 


1999-04-10 


WFPC2 


F656N 


1600 








07273_01 


1999-04-10 


WFPC2 


F814W 


40 


Circinus ULX4 


356 


25 


09379_64 


2002-12-11 


ACS-HRC 


F330W 


1200 






ACIS-S 


07273_01 


1999-04-10 


WFPC2 


F502N 


1800 








07273-01 


1999-04-10 


WFPC2 


F547M 


60 








06359_08 


1996-08-11 


WFPC2 


F606W 


600 








07273_01 


1999-04-10 


WFPC2 


F656N 


1600 








07273 _01 


1999-04-10 


WFPC2 


F814W 


40 


NGC 2403 X-l 


2014 


35 


10579_a3 


2005-10-17 


ACS-HRC 


F330W 


2912 






ACIS-S 


10579_03 


2005-10-17 


ACS-WFC 


F435W 


1248 








10579_03 


2005-10-17 


ACS-WFC 


F606W 


1248 


NGC5128 ULX1 


7797 


97 


06789_al 


1997-07-27 


WFPC2 


F555W 


480 






ACIS-I 


10260-12 


2004-8-11 


ACS-WFC 


F606W 


2370 








06789_al 


1997-07-27 


WFPC2 


F814W 


450 


CXOU J132518.3 


7797 


97 


06789_al 


1997-07-27 


WFPC2 


F555W 


480 


-430304 




ACIS-I 


10260_12 


2004-8-11 


ACS-WFC 


F606W 


2370 








06789_al 


1997-07-27 


WFPC2 


F814W 


450 


NGC 4736 XMM1 


808 


47 


10402.06 


2005-05-24 


WFPC2 


F336W 


1800 






ACIS-S 


09042_80 


2001-07-02 


WFPC2 


F450W 


460 








10402_06 


2005-05-25 


WFPC2 


F555W 


400 








09042.80 


2001-07-02 


WFPC2 


F814W 


460 


Hoi II X-l 


1564 


5 


10522.03 


2006-01-28 


ACS-WFC 


F502N 


1650 






ACIS-S 


10522_03 


2006-01-28 


ACS-WFC 


F550M 


1505 








10522_03 


2006-01-28 


ACS-WFC 


F658M 


1680 








10522.03 


2006-01-28 


ACS-WFC 


F660N 


1686 








10522_03 


2006-01-28 


ACS-WFC 


F814W 


600 


M83 XMM1 


793 


51 


10579_al 


2006-02-25 


ACS-HRC 


F330W 


2568 






ACIS-S 


10579_11 


2006-02-25 


ACS-WFC 


F435W 


1000 








10579-11 


2006-02-25 


ACS-WFC 


F606W 


1000 


M83 XMM2 


793 


51 


09774_af 


2004-07-28 


WFPC2 


F336W 


2400 






ACIS-S 


09774_0f 


2004-07-28 


ACS-WFC 


F435W 


680 








09774_0f 


2004-07-28 


ACS-WFC 


F555W 


680 








09774 _0f 


2004-07-28 


ACS-WFC 


F814W 


430 


NGC 5204 X-l 


3943 


5 


09370-01 


2002-10-29 


ACS-HRC 


F220W 


2720 






ACIS-S 


09370_01 


2002-10-28 


ACS-HRC 


F435W 


2600 








08601 _39 


2001-05-28 


WFPC2 


F606W 


600 








08601 _39 


2001-05-28 


WFPC2 


F814W 


600 


NGC 5408 X-l 


4557 


5 


08601-41 


2000-07-04 


WFPC2 


F606W 


600 






ACIS-S 


08601 Al 


2000-07-04 


WFPC2 


F814W 


600 


NGC 3034 ULX3 


8505 


83 


10776_24 


2006-03-27 


ACS-WFC 


F435W 


450 






HRC-S 


10776_24 


2006-03-27 


ACS-WFC 


F555W 


350 








10776-24 


2006-03-27 


ACS-WFC 


F658M 


1100 








10776_24 


2006-03-27 


ACS-WFC 


F814W 


175 
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Source Chandra Expos. (ks) HLA Date of Instrument Filter Expos. (s) 

ObsID ObsID Instr. HST obs 



NGC 3034 ULX4 


2933 


18 


10776_24 


2006-03-27 


ACS 


-WFC 


F435W 


450 






ACIS-S 


10776_24 


2006-03-27 


ACS 


-WFC 


F555W 


340 








10776_24 


2006-03-27 


ACS 


-WFC 


F658M 


1100 








10776_24 


2006-03-27 


ACS- 


-WFC 


F814W 


175 


NGC 3034 ULX5 


8505 


83 


10776_24 


2006-03-27 


ACS 


-WFC 


F435W 


450 






HRC-S 


10776_24 


2006-03-27 


ACS 


-WFC 


F555W 


340 








10776_24 


2006-03-27 


ACS 


-WFC 


F658M 


1100 








10776_24 


2006-03-27 


ACS 


-WFC 


F814W 


175 


NGC 3034 ULX6 


379 


9 


10776_24 


2006-03-27 


ACS 


-WFC 


F435W 


450 






ACIS-I 


10776_24 


2006-03-27 


ACS 


-WFC 


F555W 


340 








10776_24 


2006-03-27 


ACS 


-WFC 


F658M 


1100 








10776_24 


2006-03-27 


ACS 


-WFC 


F814W 


175 


CXOU J095550.6 


10542 


118 


10776_24 


2006-03-27 


ACS 


-WFC 


F435W 


450 


+694044 




ACIS-S 


10776_24 


2006-03-27 


ACS 


-WFC 


F555W 


340 








10776_24 


2006-03-27 


ACS 


-WFC 


F658M 


1100 








10776_24 


2006-03-27 


ACS 


-WFC 


F814W 


175 
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Table 3 

Details of observations used to locate possible optical counterparts to our ULX 

sample. 



Source 


No. of Matches 


ULX Position 


Filter^ 


Positional error 




matches' 2 used?^ 


RA 


DEC 




RA e (*) 




NGC 598 ULX1 


N 


01 33 50.90 


+30 39 36.76 


F170W 


0.79 


0.81 










F336W 


0.81 


0.88 










F439W 


0.91 


0.91 










F555W 


0.78 


0.78 










ThQ 1 AWT 


0.68 


0.71 


NGC 55 ULX1 


1 N( J > 


00:15:28.90 


-39:13:18.77 


F606W 


0.71 


0.71 










F814W 


0.71 


0.66 




y( j ) 


00-1 5-9R QO 


— D 7.1 J.IO. / O 


F606W 


0.45 


0.45 










F814W 


0.46 


0.37 


NGC 4190 X-l 


N 


12 13 45.27 


+36 37 54.66 


F300W 


0.55 


0.55 










F450W 


0.58 


0.58 










F606W 


0.56 


0.56 










F814W 


0.57 


0.57 


NGC 253 ULX1 


N 


00 47 34.00 


-25 16 36.35 


F336W 


0.65 


0.71 










F475W 


0.61 


0.63 










F555W 


0.77 


1.09 










F606W 


0.61 


0.62 










F814W 


0.60 


0.62 


NGC 253 ULX2 


2 N<2) 


00 47 32.97 


-25 17 48.92 


F475W 


0.61 


0.63 










F606W 


0.61 


0.62 










F814W 


0.60 


0.62 


NGC 253 ULX3 


2 N<2) 


00 47 33.44 


-25 17 21.99 


F336W 


0.63 


0.71 










F475W 


0.61 


0.63 










F502N 


0.71 


0.87 










F555W 


0.59 


1.09 










F606W 


0.61 


0.62 










F814W 


0.60 


0.62 


NGC 253 XMM6 


2 N^) 


00 47 42.77 


-25 15 02.15 


F475W 


0.61 


0.63 










F606W 


0.59 


0.62 










F814W 


0.57 


0.60 


M81 X-6 


2 Y 


09 55 32.94 


+69 00 33.66 


F336W 


0.42 


0.44 










F435W 


0.41 


0.44 










F555W 


0.42 


0.44 










F606W 


0.42 


0.44 










F814W 


0.42 


0.44 


Hoi IX X-l 


2 Y 


09:57:53.28 


+69:03:48.31 


F330W 


0.39 


0.38 










F435W 


0.39 


0.37 










F555W 


0.25 


0.38 










F814W 


0.54 


0.37 


NGC 4395 UI X1 


1 N< 5 ) 


12 26 01.44 


+33 31 30.99 


F336W 


0.77 


0.69 










F435W 


0.50 


0.58 










F555W 


yj.jyj 


U.JO 










F814W 


0.50 


0.58 


NGC 1313 X-l 


3 Y 


03 18 20.02 


-66 29 10.85 


F330W 


0.47 


0.46 










F435W 


0.46 


0.46 










F555W 


0.48 


0.40 










F606W 


0.28 


0.25 










F814W 


0.52 


0.38 


NGC 1313 X-2 


4 Y 


03 18 22.24 


-66 36 03.61 


F330W 


0.95 


0.76 










F435W 


0.32 


0.36 










F555W 


0.42 


0.32 










F814W 


0.44 


0.36 


IC 342 X-l 


4 Y 


03:45:55.63 


+68:04:55.42 


F330W 


0.58 


0.62 










F435W 


0.24 


0.36 










F555W 


0.21 


0.35 










F606W 


0.21 


0.36 










F658M 


0.23 


0.36 










F814W 


0.22 


0.35 
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Source 


No. of Matches 


ULX Position 


Filter^ 


Positional error 




matches" used?^ 


RA 


DEC 






Dec y ' 


IC 342 X-2 


2 N<^) 


Af\ 1 S 7^ 
UJ ^tO 1j. / J 


too 1 1 1Z.OJ 


F330W 


U.OJ 


n M 










F435W 


0.64 


0.69 










F606W 


0.65 


0.70 




9 vU) 


03 46 15.87 


+68 11 12.92 


rjjUW 


0.25 


0.60 










F435W 


0.25 


0.60 










F606W 


0.27 


0.61 


IC 342 ULX2 


4 N' 5 ) 


03 46 48.52 


+68 05 46.83 


F555W 


0.67 


0.67 










F606W 


0.69 


0.87 










F656N 


0.62 


0.67 










F675W 


0.66 


0.71 










F814W 


0.61 


0.67 


TC 34? X-6 


4 Y 


(W-Af^-'M zLI 
Uj.^+O. J / .^tl 


tOo.UO. 1 o.oO 


F555W 


U.J 1 


u.oz 










F606W 


U. JO 


u.oo 










F814W 


n aa 


U.OJ 


Circinus ULX1 


4 N^) 


14 13 12.22 


-65 20 13.85 


F502N 


0.55 


0.57 










F547M 


0.63 


0.64 










F606W 


0.56 


0.57 










F656N 


U. Jy 


u.oz 










F814W 


U.J J 


U.JO 


Circinus ULX3 


4 N^) 


14 13 10.26 


-65 20 17.97 


F330W 


0.55 


0.57 










F502N 


0.55 


0.57 










F547M 


0.63 


0.68 










F606W 


0.56 


0.57 










F656N 


U. Jy 


n 69 
u.oz 










F814W 


U.J J 


U.JO 


Circinus ULX4 


4 


1 A 1 ^ 1 ^ 


on 99 A*> 
— OJ ZU ZZ.Hj 


F330W 


U.J J 


n S7 

U.J / 










F502N 


U.J J 


n S7 

U.J / 










F547M 


0.63 


0.68 










F606W 


0.56 


0.57 










F656N 


0.59 


0.62 










F814W 


0.55 


0.56 


NGC 2403 X-l 


i n( ? ) 


07 36 25.57 


+65 35 39.88 


F330W 


0.65 


0.70 










F435W 


0.65 


0.70 










F606W 


0.64 


0.64 




y( 7 ) 


07:36:25.52 


+65:35:40.01 


F330W 


0.51 


0.57 










F435W 


0.51 


0.57 










F606W 


0.50 


0.50 


NGC5128 ULX1 




13 25 19.83 


-43 03 16.20 


F555W 


0.81 


0.84 










F606W 


0.72 


0.73 










F814W 


0.82 


0.82 




y( § ) 


13 25 19.88 


-43 03 16.25 


F555W 


0.62 


0.66 










F606W 


0.50 


0.57 










F814W 


0.60 


0.60 


CXOUJ1 325 18.3 




13 25 18.31 


-43 03 04.63 


F555W 


0.70 


0.81 


-430304 








F606W 


0.67 


0.69 










F814W 


u. / 1 


n 87 

U.o / 




y(^) 


13 25 18.24 


-43 03 04.50 


F555W 


0.58 


0.63 










F606W 


0.45 


0.52 










F814W 


0.58 


0.58 


NGC 4736 XMM1 


1 N^) 


12 50 50.33 


+41 07 12.19 


F336W 


0.85 


0.77 










F450W 


0.81 


0.82 










F555W 


0.83 


0.77 










F814W 


0.74 


0.69 


Hoi II X-l 


N 


08 19 29.00 


+70 42 19.08 


F502N 


0.58 


0.54 










F550M 


0.54 


0.53 










F658M 


0.54 


0.53 










F660N 


0.63 


0.53 










F814W 


0.54 


0.53 


M83 XMM1 


N 


13 37 19.80 


-29 53 48.80 


F330W 


0.55 


0.54 










F435W 


0.55 


0.54 










F606W 


0.56 


0.56 
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Source 


No. of 
matches a 


Matches 
used? 6 


ULX Position 
RA DEC 


Filter^ 


Positional error 

RA e W Dec e (^) 


M83 XMM2 





N 


13 36 59.45 


-29 49 59.21 


F336W 
F435W 
F555W 
F814W 


0.76 
0.62 
0.59 
0.58 


0.77 
0.67 
0.61 
0.67 


NGC 5204 X-l 





N 


13 29 38.61 


+58 25 05.55 


F220W 
F435W 
F606W 
F814W 


1.04 
1.07 
0.91 
0.92 


1.30 
1.09 
0.94 
0.87 


NGC 5408 X-l 





N 


14 03 19.61 


-41 22 58.65 


F606W 
F814W 


0.54 
0.54 


0.60 
0.60 


NGC 3034 ULX3 





N 


09 55 51.33 


+69 40 43.65 


F435W 
F555W 
F658M 
F814W 


0.67 
0.64 
0.60 
0.61 


0.60 
0.59 
0.59 
0.61 


NGC 3034 ULX4 





N 


09 55 51.02 


+69 40 45.02 


F435W 
F555W 
F658M 
F814W 


0.67 
0.64 
0.60 
0.61 


0.60 
0.60 
0.59 
0.60 


NGC 3034 ULX5 





N 


09 55 50.17 


+69 40 46.47 


F435W 
F555W 
F658M 
F814W 


0.67 
0.64 
0.60 
0.61 


0.60 
0.60 
0.59 
0.61 


NGC 3034 ULX6 





N 


09 55 47.46 


+69 40 36.28 


F435W 
F555W 
F658M 
F814W 


0.75 
0.72 
0.69 
0.69 


0.69 
0.69 
0.68 
0.69 


CXOU J095550.6 
+694044 





N 


05 55 50.65 


+69 40 43.81 


F435W 
F555W 


067 
0.64 


0.69 
0.60 



F658M 0.60 0.59 
F814W 0.61 0.61 



Note. — a Number of cross-correlated sources found within 6' of the ULX. * denotes position improved via private communication with 
authors of previous paper, see note (8) & (9) for more details. b Have cross-correlated sources been used to obtain the ULX's position and 
the corresponding positional error region; yes (Y), or no (N). c Position of ULX using WAVDETECT , in some cases corrected via relative 
astrometry. d HST filter band. e Combined (relative & astrometric) error of the (/) right ascension and (//) declination of the ULX in the 
corresponding HST filter band. Notes pertaining to individual sources: (i) Field contains only one cross-matched counterpart, so caution must 
be taken. Positional error is calculated using Chandra's absolute positional accuracy in the first instance (N) and using relative astrometry in 
the second (Y). ^Cross-matched sources are in a crowded galaxy centre. Derived errors are as large as they would be without correction and 
so relative astrometry is not applied. ^Only available cross-correlated source is the nucleus. Astrometry is not corrected, due to the increased 
uncertainty in matching galaxy centres. ^ Source was not present in the deepest observation of this field, so position was matched to that of 
two other ULXs in the field. ^Chandra data is too shallow to allow accurate relative astrometric corrections. ^Jackknife errors strongly 
suggest that the 2MASS matches in the Chandra fields are incorrect, and they also produce similarly large (0.5") errors. Therefore, we opt 
not to incorporate these into our error calculations, and use the standard Chandra uncertainty instead. ^Only one cross-correlated source 
available so positions & errors are calculated using both techniques. ^ Higher accuracy achieved in Woodley et al. (2008), with position 
checked and refined via private communication with Sivakoff (2012). (9) Higher accuracy achieved in Sivakoff et al. (2008), with position 
checked and refined via private communication with Sivakoff (2012). (i0) Only possible match is from confused region around the centre of 
the galaxy. 
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Table 4 

Potential optical counterparts for each ULX, with their corrected magnitudes. 



Source 


Filter 


C/P 

1 1 J 


Previous 

TF>k 

1 1 j 


Observed 
Vega mag c 


Galaxy corrected 
Vega mag^ 


Estirnated intrinsic 
Vega niag e 


NGC 598 ULX1 


F170W 







- 


- 


- 




F336W 














F439W 














F555W 














F814W 






- 


- 


- 


NGC 55 ULX1 


F606W 


1 


N 


23.2 ± 0.3 


23.2 ± 0.3 


22 ± 2 




F814W 






23.2 ± 0.5 


23.1 ± 0.5 


22 ± 2 




F606W 


2 


N 


24.4 ± 0.6 


24.3 ± 0.6 


23 ± 2 




F814W 






23.1 ± 0.5 


23.1 ± 0.6 


22 ± 2 


NGC 4190 X-l 


F300W 


1 


N 


23.2 ± 0.3 


93 o + 3 


22 ± 4 




F450W 






24 ± 2 


24 ± 2 


23 ± 4 




F606W 






24 ± 2 


24 ± 2 


24 ± 4 




F814W 






25 ± 3 


25 ± 3 


25± 5 


NGC 253 ULX1 


F336W 


1 


N 










F475W 






26 ± 2 


26 ± 2 


20 ± 12 




F555W 














F606W 














F814W 














F336W 


2 


N 










F475W 






25 ± 2 


25 ± 2 


20 ± 12 




F555W 














F606W 






23.6 ± 0.7 


23.5 ± 0.7 


20 ± 12 




F814W 














F336W 


3 


N 










F475W 






26 ± 3 


26 ± 3 


21 ± 13 




F555W 














F606W 






25.0 ± 0.7 


25.0 ± 0.7 


21 ± 13 




F814W 






- 


- 


- 




F336W 


4 


N 










F475W 






26 ± 3 


26 ± 3 


20 ± 13 




F555W 














F606W 






25 ± 2 


25 ± 3 


21 ± 13 




F814W 












NGC 253 ULX2 


F475W 


1 


N 


21.8 ±0.3 


21.7 ±0.3 


21 ±2 




F606W 






20.9 ± 0.2 


20.8 ± 0.2 


20 ± 2 




F814W 






19.7 ± 0.2 


19.6 ± 0.2 


19 ± 2 


NGC 253 TJT X3 


F336W 


1 


N 










F475W 






22.9 ± 0.7 


22 8 + 7 


21 ± 12 




F502N 














F555W 






23 ± 2 


23 ± 2 


21 ± 12 




F606W 






22 3 + 6 


22.2 ± 0.6 


21 ± 12 




F814W 






21.7 ± 1.0 


21.6 ± 1.0 


21 ± 12 


NGC 253 XMM6 


F475W 


1 


N 


26 ± 2 


26 ± 2 


24 ± 2 




F606W 






25 ± 1 


25 ± 1 


24 ± 2 




F814W 






23 ± 4 


23 ± 4 


22 ± 5 




F475W 


2 


N 










F606W 






25.4 ± 0.8 


25.3 ± 0.8 


24 ±2 




F814W 






23.6 ± 0.9 


23.6 ± 0.9 


23 ±2 




F475W 


3 


N 










F606W 






_ 


_ 


_ 




F814W 






25 ±3 


23 ±3 


24 ±4 




F475W 


4 


N 










F606W 






- 


- 


- 




F814W 






24 ± 2 


24 ± 2 


23 ± 2 




F475W 


5 


N 










F606W 






- 


- 


- 




F814W 






24 ± 1 


24 ± 1 


23 ±2 


M81 X-6 


F336W 


1 


Y (2,S) 


22 ±2 


22 ±2 


21 ±3 




F435W 






24.1 ± 0.7 


23.7 ± 0.6 


23 ±2 




F555W 






24 ±2 


24 ±2 


23 ±3 




F606W 






24.1 ± 0.6 


23.8 ± 0.6 


23 ±2 




F814W 






23 ±2 


23 ±2 


22 ±3 


Hoi IX X-l 


F330W 


1 


Y (5) 


21.0 ±0.5 


20.6 ± 0.5 


20 ± 1 




F435W 






22.6 ± 0.3 


22.3 ± 0.3 


22 ± 1 




F555W 






22.8 ± 0.4 


22.5 ± 0.3 


22 ± 1 
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Source 


Filter 


C/P 


Previous 


Observed 


G3.l3.xy corrected 


Estirri3ted intrinsic 






1 Y J 


1 1 j 


Veg3. ni3g c 


Veg3 m3g^ 


Veg3 ni3g e 




F814W 






22.3 zb 0.3 


22.2 zb 0.3 


22 zb 1 




F330W 


2 


N 


24 ±2 


24zb2 


22zb2 




F435W 














F555W 






_ 


_ 


_ 




F814W 






- 


- 


- 




F330W 


3 


N 










F435W 






26 ±2 


26zb2 


25 zb2 




F555W 






21 ±2 


26zb2 


26zb3 




F814W 






26 ±2 


26zb2 


26zb2 


NGC 4395 ULX1 


F336W 


1 


N 


22 ± 1 


22 zb 1 


21 zb3 




F435W 






CG 


CG 


CG 




F555W 






CG 


CG 


CG 




F814W 






CG 


CG 


CG 


NGC 1313 X-l 


F330W 


1 




23 ± 1 


22 zb 1 


21 zb 1 




F435W 






24.1 ± 1.0 


23.6 zb 1.0 


23zbl 




F555W 






24 ±2 


24 zb 1 


23 zb2 




F606W 






24.0 zb 0.8 


23.7 zb 0.8 


23.0 zb 0.9 




F814W 






24 ± 1.0 


24.0 zb 1.0 


23zbl 


NGC 1313 X-2 


F330W 


1 


y(5,5) 


22.0 zb 0.8 


21.6 zb 0.8 


20zb2 




F435W 






23.5 zb 0.4 


23.2 zb 0.4 


22zb2 




F555W 






23.6 ± 0.5 


23.4 =b 0.5 


22zb2 




F814W 






23.6 zb 0.6 


23.5 zb 0.6 


23 zb2 




F330W 


2 


N 


- 


- 


- 




F435W 






- 


- 


- 




F555W 






24.3 zb 0.7 


24.1 zb0.7 


23 zb2 




F814W 






- 


- 


- 


IC 342 X-l 


F330W 


1 


Y( e > 7 ) 


_ 


_ 


_ 




F435W 






25.2 zb 1.0 


23.1 zb 1.0 


22zb2 




F555W 






24.1 ±0.6 


22.3 zb 0.6 


20 zb 1 




F606W 






23.6 zb 0.4 


21.9 zb 0.4 


20 zb 1 




F658M 






23. ±2 


22zb2 


20zb2 




F814W 






22.2 ± 0.3 


21.1 zb0.3 


20 zb 1 




F330W 


2 


N 


- 


- 


- 




F435W 






27 ±2 


25zb2 


22zb2 




F555W 






26 ±2 


24zb2 


22zb2 




F606W 






26 ±2 


24zb2 


22zb2 




F658M 






- 


- 


- 




F814W 






25 zb 1 


24 zb 1 


22zb2 


IC 342 X-2 


F330W 


1 


\ \ 1 ) 










F435W 






27 ± 1 


25 zb 1 


25zb3 




F606W 














F330W 


2 


N 


_ 


_ 


_ 




F435W 






28 ±2 


25 zb2 


25zb3 




F606W 














F330W 


3 


N 










F435W 






21 ±2 


25 zb2 


25 zb2 




F606W 












IC 342 ULX2 


F555W 





N 


_ 


_ 


_ 




F606W 






_ 


_ 


_ 




F656N 






_ 


_ 


_ 




F675W 














F814W 






_ 


_ 


_ 


IC 342 X-6 


F555W 


1 


N 


26 zb 4 


24 zb 3 


24 zb 9 




F606W 














F814W 






26 ±2 


24. zb2 


24=b6 




F555W 


2 


N 


- 


- 


- 




F606W 






- 


- 


- 




F814W 






25 zb 1 


24 zb 1 


25zb6 




F555W 


3 


N 










F606W 














F814W 






25.7 zb 0.9 


24.6 zb 0.9 


25zb6 


Circinus ULX1 


F502N 


1 


y(«) 










F547M 














F606W 






24=b6 


20zb5 


20zb5 




F656N 














F814W 
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Table 4 — Continued 



Source 


Filter 


C/P 


Previous 


Observed 


Galaxy corrected 


Estimated intrinsic 






ID a 


lD b 


Vega mag c 


Vega mag d 


Vega mag e 


Circinus ULX3 


F330W 





N 


_ 


_ 


_ 




F502N 






_ 


_ 


_ 




F547M 






_ 


_ 


_ 




F606W 






>22 


- 


- 




F656N 














F814W 






_ 


_ 


_ 


Circinus ULX4 


F330W 





N 


- 


- 


- 




r juzin 














F547M 






- 


- 


- 




F606W 






> 22 


- 


- 




F656N 






— 


— 


— 




F814W 












NGC 2403 X-l 


F330W 


1 


y( 7 ) 


24 ±2 


23 ±2 


21 ±3 




F435W 






25 ±1 


25 ±1 


23 ±3 




rOUOW 






OA fx A- H f 
Z"4-.0 Zt U.O 


OA ^ 4- n f\ 
Z"4-.J zt U.O 


Zj Zt J 


NGC5128 ULX1 


F555W 





y(*) 


>24 


_ 







F606W 








_ 







F814W 






- 


- 


- 


CXOU J132518.3 


F555W 





N 


>24 


_ 


_ 


-430304 


F606W 






- 


- 


- 




F814W 






- 


- 


- 


NGC 4736 XMM1 


F336W 





N 


- 


- 


- 


















F555W 






>24 


_ 


_ 




F814W 






- 


- 


- 


Hoi II X-l 


F502N 


1 


-\r( 10) 










F550M 






21.6 ±0.3 


21.5 ±0.3 


21 ±2 




F658M 














F660N 






21 ±2 


21 ± 1 


21 ±2 




F814W 






21.6 ±0.2 


21.5 ±0.2 


21 ±2 


M83 XMM1 


F330W 


1 


N 


_ 


_ 


_ 




F435W 






26 ±3 


26 ±3 


25 ±9 




F606W 






26 ±2 


25 ±2 


25 ±8 


M83 XMM2 


F336W 





N 


>21 


- 


- 




F435W 






CG 


CG 


CG 




F555W 






CG 


CG 


CG 




F814W 






CG 


CG 


CG 


NGC 5204 X-l 


F220W 


1 


y( jj ) 


20.0 ± 0.5 


19.8 ±0.5 


19 ±3 




F435W 






22.40 ±0.10 


22.37 ±0.10 


22 ±3 




F606W 






22.3 ± 0.8 


22.3 ± 0.8 


22 ± 3 




F814W 






23 ± 1 


23 ±1 


23 ±3 




F220W 


2 




19.6 ± 0.5 


19.5 ± 0.5 


19 ±3 




F435W 






20.8 ± 0.2 


20.8 ±0.1 


21 ±3 




F606W 






20.0 ±0.1 


20.0 ±0.1 


20 ±3 




ThQ 1 AWT 






1 Q f -U n 1 
ly.D zt U.Z 


ly.O zt U.Z 


in j. o 
ly zt J 


NGC 5408 X-l 


F606W 


1 


Y (l3) 


22.4 ± 0.6 


22.2 ± 0.6 


22 ±2 




F814W 






23 ±3 


23 ±3 


23 ±4 


NGC 3034 ULX3 


F435W 





N 










F555W 






^20 








F658M 








_ 


_ 




F814W 






_ 


_ 


_ 


NGC 3034 ULX4 


F435W 





N 










F555W 






^23 








F658M 














F814W 












NGC 3034 ULX5 


F435W 


1 


N 


24 ±2 


23 ±2 


17 ±2 




F555W 






22.3 ± 0.8 


21.8 ±0.8 


17.0 ±0.8 




F658M 














F814W 













44 Gladstone et al. 



Table 4 — Continued 



Source 


Filter 


C/P 


Previous 


Observed 


G3.l3.xy corrected 


Estim3ted intrinsic 






1 Y J 


1 u 


Veg3. ni3g c 


Veg3 ni3g^ 


Veg3 ni3g e 




F435W 


2 


N 


22.6 ± 0.7 


21.9 ± 0.7 


15.7 ± 0.3 




F555W 














F658M 






_ 


_ 


_ 




F814W 






_ 


_ 


_ 


NGC 3034 ULX6 


F435W 


1 


N 










F555W 














F658M 














F814W 






18.46 ± 0.07 


18.16 ±0.07 


16* 


CXOU J095550.6 


F435W 





N 








+694044 


F555W 






^23 








F658M 














F814W 













Note. — "Candidate counterpart identification number, is listed if no counterpart is available b Previous identification as the candidate 
optical counterpart to this source, references as follows: (i) Dubus et al. (2004), (s) Liu et al. (2002), (s) Ramsey et al. (2006), (5) Grise et al. 
(2005), U) Yang et al. (2011), (5) Liu et al. (2007), (e) Feng & Kaaret (2008), (7) Roberts et al. (2008), (5) Weisskopf et al. (2004), (9) Ghosh et 
al. (2006), (i0) Kaaret et al. (2004), (ii) Liu et al. (2004), (iS) Goad et al. (2002), (i5) Lang et al. (2007). c Aperture corrected observed optical 
magnitude of each potential counterpart. CG denotes cases where error region resides in a chip gap. ~ my limits are provided in those cases 
where no counterpart was observed. M83 XMM2 is located in the chip gap in the ~ V band, so a limit is obtained for the only band where 
this is not the case - F336W. d Galactic extinction corrected Vega magnitudes. 6 Estimated intrinsic Vega-mag for each potential counterpart, 
calculated using published values of Nh (provided in TablefTJ, assuming the Milky Way value of R v . 
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Table 5 

Typing of potential ULX counterparts. 



Source 


C/P 


Type b 


Type 


Distance 




offset^ 


Am 5 55 9 


Types from 


M 55 5* 


Types from 


Previous 




1 1 J 




A/fx r C 

IVly 


modulus'^ 








colour^ 1 




A/7^^ 3 


1 VJ 


NGC 598 ULX1 


nucleus 


— 


— 


24.771 


— 


— 


— 


— 


— 


— 


nucleus 


NGC 55 ULX1 


1 


Al V 


+0.7* (*) 


26.203 


26.898 


23.179 


-3.719 


O, B, A, 


-3.0 


late B V, 


- 


















early F 




IV, III, A II 






2 


MOV 


+9.0^) 


26.203 


35.208 


24.704 


-10.499 


late K, M 


-1.499 


B V, IV, 


- 


















early F 




& late III, 




















early F 




F-K II 




NGC 4190 X-l 


1 


B2 V 


-2.641 («) 


27.236 


24.595 


24.114 


0.481 


O, B, A, 


-3.1 


mid B III, 


- 


















F, G, early 




late A or 




















to mid F 




later II 




NGC 253 ULX1 


1 






27.386 














- 




2 


M2 III 


-0.6< J > 


27.386 


29.930 


24.099 


-2.831 


all 


-3.3 


B V, late 
























B III, A II 






3 


K0 III 


+0.5^) 


27.386 


27.886 


25.234 


-2.652 


all 


-2.152 


LateB 


- 






















all B III, 
























all II 






4 


G8 III 


+0.6^) 


27.386 


28.006 


25.073 


-2.933 


all 


-2.313 


LateB 


- 






















all B III, 
























all II 




NGC 253 ULX2 


1 


MOV 


+9.0^) 


27.386 


36.478 


21.165 


-15.313 


M, late K 


-6.2 


all lb 


- 


NGC 253 ULX3 


1 


K8 V 


+6.7* <*> 


27.386 


34.126 


22.518 


-11.608 


Late B, A, 


-4.868 


OB V & III, 


- 


















F G K 




A II 




















early M 




A-M lb 




NGC 253 XMM6 


1 


MOV 


+9.0( J ) 


27.386 


36.478 


25.601 


-10.877 


all 


-1.8 


Late B V 


- 






















& IV, F-M II 






2 


MO III 


-0.2( J ) 


27.386 


27.304 


25.782 


-1.522 


M 


-1.4 


Late B V 
























& IV, F-M II 






3 


- 


- 


27.386 


- 


- 


- 


- 


- 


- 


- 




4 


- 


- 


27.386 


- 


- 


- 


- 


- 


- 


- 




5 


- 


- 


27.386 


- 


- 


- 


- 


- 


- 


- 


IVlol 


i 
i 


R9 TTT 
£>Z 111 


-Z.OJ> v y 


L 1 .OJ / 


LJ XjL 1 


Zj. /yj 


-Y.LjL 


vJ. 15, r\, 


-j.DoZ 


R V TV 
D V, IV, 


Wo v v ^ ' ' 


















F, 




Late B II, 




















early G^ 




A II 




Hol IX X-l 


1 


B4 V 


-1.22(«) 


27.670 


26.429 


22.253 


-4.176 


B 


-5.4 


O, early 


OB( 5 ' e ) 






















B V BII 
























all lb ' 






2 






27.670 














— 




3 


B6 V 


-0.8l( s ) 


27.670 


zu.ouu 


26.123 


-0.737 


O B A 


-1.574 


T atp R V 




















F,'mi'dG 




& II, F-K II 




M/" 1 /IIOC TTT V1 

NCjC 439j ULX1 


1 






T7 TOO 
















NGC 1313 X-l 


1 


B2 III 


-2.63^ j 


27.841 


25.180 


23.707 


-1.473 


O, B 


-4.1 


Late O V, 


K5-M0 IL ' > 


















F, mid G 




IV, III, B, lb 




NGC 1313 X-2 


1 


B2 III 


-2.63^ j 


27.841 


25.180 


23.359 


-1.821 


O, early B 


-4.5 


OB, all lb 






2 






27.841 
















IC 342 X-l 


1 


K0 IV 


+3.2( J ) 


27.955 


31.221 


22.227 


-8.994 


K, G 


-5.7 


OB, all I 


F8-G0 IbW, 
























F0-F5 l( i0 ) 




2 


F6 V 


+3 .4* (J) 


27.955 


31.378 


24.045 


-7.333 


O, B, A, 


-3.9 


OB & A II 




















F, G, K, 




Fib 




















early M 




or later 




IC 342 X-2 


1 






27.955 














y(10) 




2 






27.955 


















3 






27.955 
















IC 342 ULX2 


nucleus 






27.955 
















IC 342 X-6 


1 


B6V 


-0.8l(^) 


27.955 


27.145 


24.407 


-2.738 


O, B, A, 


-3.548 


B V, IV, 




















F, G, K, 




III, A II, 
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Table 5 — Continued 



Source 


C/P 
ID a 


Type b 


Type 


Distance 
modulus d 


m 5 55 e 


offset^ 


Am 5 55 9 


Types from 
colour^ 


M555* 


Types from 

M 55 5 J ' 


Previous 
ID fc 




2 
3 


; 




27.955 
27.955 


; 


; 


: 


early M 


; 


KIb 
or later 


; 


Circinus ULX1 


1 


_ 


_ 


28.010 


_ 


_ 


_ 


_ 


_ 


_ 


K5 or 

later(^) 


NGC 2403 X-l 


1 


B8Ia 


-5.67( 2 ) 


28.116 


22.444 


24.538 


+2.095 


O, B, A, 
F, 

early K 


-3.6 


B V, IV, 
III, A II, 

KIb 
or later 


OB giant/ 
supergiant( i0 ) 


NGC 5128 ULX1 









28.121 














OB (^) 


Hoi II X-l 


1 


A3 III 


-0.2*( J ) 


28.266 


28.756 


21.490 


-7.266 


B, A 


-6.8 


all la 


04 V/ 
B3 Ib^ ; 


M83 XMM1 


1 


A3 V 


+2.0* (*) 


28.360 


30.358 


25.486 


-4.872 


all 


-2.9 


Late B, 

All 
or later 




NGC 5204 X-l 


1 


05 V 


-5.00( s ) 


28.406 


23.370 


22.743 


-0.627 


O 


-5.3 


0,BII 


05 V, 
07III or 

BOIb(^) 




2 


F8 V 


+3.4* (*) 


28.406 


31.840 


20.184 


-11.656 


none 


-8.2 


Late 
A la, 
FG la or 
early 
Mia 


star cluster^ 5 ) 
05 V + 
cluster^ ) 


NGC 5408 X-l 


1 


08f 


-4.83* (2) 


28.406 


23.539 


22.198 


-1.341 


O, B, A, 
F, G, K, 
early M 


-6.2 


o, 

BII&I 


B/Al(^) 


NGC 3034 ULX5 


1 

2 


M2 III 


-0.6(^) 


28.580 
28.580 


28.108 


21.610 


-6.498 


B, A, F, G, 
K,M 


6.9 


all la 





NGC 3034 ULX6 1 - - 28.580 



Note. — a Candidate counterpart ID taken from Table [4] We remove those ULXs for which there are no candidate counterparts detected 
from the list. The only exception to this is NGC 5128 ULX1, as a candidate counterpart has been published elsewhere. b Stellar type derived 
from x 2 fitting of HST band stellar templates, created using SYNPHOT. c Absolute magnitude of identified stellar type in V band (taken from 
similar type if particular My is unavailable, sources are indicated by *). ^Calculated distance modulus using values in Table[l] e Derived Vesa 
magnitude in F555W band, assuming stellar type classification is correct, ^offset required to fit observed data to stellar types, this can also be 
considered to be the apparent magnitude in the F555W band (~ 771555). This is found by % 2 fitting. 9 Value A771555 = offset — 771555- h All 
possible types allowed from SED fitting, within errors. 1 Galactic extinction corrected absolute magnitude in HST band F555W, calculated 
using offset value. J A11 possible types allowed from absolute magnitude fitting, within errors. k Stellar types collated from work of other 
authors, where Y refers to identification without typing. Figures shown in brackets relate to following references: (^Zombeck (1990), 
(s) Wegner (2006), (5) Dubus et al. (2004), ( ^Liu et al. (2002), (5) Ramsey et al. (2006), (6) Grise et al. (2005), (7) Yang et al. (201 1), (5) Liu et 
al. (2007), (9) Feng & Kaaret (2008), (i0) Roberts et al. (2008), (ii) Weisskopf et al. (2004), (^Ghosh et al. (2006), (i5) Kaaret et al. (2004), 
(^)Liu et al. (2004), (i5) Goad et al. (2002), (^Lang et al. (2007). 
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Table 6 

Models, flux ratios and X-ray luminosities for each ULX that has a counterpart 
detected in multiple HST bands. 



Source 


Model a 


Flux ratio 6 


If 












(10 39 er^ 


o ) 






Galactic 


Intrinsic 


Galactic 


Intrinsic 


NGC 55 ULX1 


DKBBFTH^) 


2.60 


0.873 


1.1 


2.0 


NGC 4190 X-l 


diskpbb( s ) 


6.80 


3.17 


4.0 


4.7 


NGC 253 ULX1 


PILEUP x PEGPWRLW ( * ) 


104 


1.29 


0.64 


1.6 


NGC 253 ULX2 


diskbb( 5 ) 


13.1 


6.62 


1.5 


1.9 


NGC 253 ULX3 


POWERLAW^*) 


33.1 


7.10 


0.19 


0.23 


NGC 253 XMM6 


PILEUP x (DISKBB+COMPTT)^ 


15.3 


2.85 


0.93 


1.5 


M81 X-6 


DKBBFTH^) 


9.33 


3.93 


2.3 


2.7 


Hoi IX X-l 


DKBBFTH^) 


8.08 


4.25 


7.7 


8.7 


NGC 1313 X-l 


DKBBFTH^) 


6.96 


2.33 


3.7 


4.8 


NGC 1313 X-2 


DKBBFTH^) 


12.4 


4.85 


4.8 


5.6 


IC 342 X-l 


DKBBFTH^) 


27.4 


3.76 


3.0 


4.2 


IC 342 X-6 


DISKBB^) 


6.34 


6.34 


0.16 


0.16 


NGC 2403 X-l 


DKBBFTH^) 


8.38 


3.07 


2.5 


3.2 


Hoi II X-l 


DKBBFTH^) 


2.05 


1.27 


15.4 


19.1 


M83 XMM1 


diskpn+eqpair( 5 ) 


0.359 


0.236 


1.0 


1.8 


NGC 5204 X-l 


DKBBFTH^) 


2.11 


1.61 


2.0 


2.2 


NGC 5408 X-l 


DKBBFTH^) 


0.685 


0.552 


4.6 


5.4 


NGC 3034 ULX5 


POWERLAW^) 


125 


3.48 


27 


42 



Note. — "Published model taken from references denoted by superscript number, although not stated, all models are absorbed. 6 Galactic 
corrected and intrinsic flux ratio (1.0 - 10.0 / 0.3 - 1.0 keV), derived within XSPEC using listed models. c Galactic absorption corrected and 
intrinsic X-ray luminosity in the 0.3 -10.0 keV bandpass. These luminosities are taken from references given below whenever the bandpass 
matches. If this was not available, we scale using model fits within XSPEC. References: ^Gladstone et al. (2009), ^Roberts et al. (in prep), 
(5) Kajavaetal. (2009), U) Maketal. (2011), (5) Stobbart et al. (2006), (e) Kaaret et al. (2006). ( * } Values not available from published results, 
so simple fitting performed on data using listed models. 
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Table 7 

la binary parameter constraints obtained for each potential ULX optical 
counterpart derived from Galactic absorption corrected magnitudes. 









M B h = 


= 10M o 


M B h = 


100M o 


M B h = 


1000M© 


Source 


C/P ID a 


cos(/) 








R* 












(Mq) 


(Rq) 


(Mq) 


(Rq) 


(Mq) 


(Rq) 


NGC 55 ULX1 


1 


0.0 


1.4 - 15.4 


3.4 - 64.1 


3.8 - 19.7 


4.1 - 59.1 


3.8 - 20.4 


4.1 - 58.9 






0.5 


< 13.9 


2.2 - 5.1 


< 8.9 


0.9 - 3.1 


A f D TT <" 

1V1 BH ^ 


700 Mr, 




2 


0.0 


1.4-1.5 


50.9-52.3 


1.4-4.6 


29.3-55.4 


1.4-4.6 


29.3-55.5 






0.5 


< 15.1 


3.6-5.4 


< 10.0 


1.4-3.3 


M B h < 


84OM 


NGC4190 X-l 


1 


0.0 


1.3-48.8 


2.0-85.9 


2.9-58.2 


3.0-439.9 


2.9-63.3 


3.0-447.6 






0.5 


< 45.9 


< 13.0 


< 43.6 


< 8.3 


< 13.4 


< 4.0 


NGC 253 ULX1 


2 


0.0 


1.4-29.0 


3.5-95.7 


2.5-32.9 


4.1-263.7 


2.6-34.4 


4.2-268.9 






0.5 


< 30.9 


2.1-12.2 


< 28.0 


0.8-6.7 


< 10.0 


< 3.3 




3 


0.0 


1.2 - 14.1 


2.1 - 84.7 


1.2 - 17.9 


2.7 - 103.2 


1.2 - 18.5 


2.7 - 103.4 






0.5 


< 12.2 


< 4.6 


< 7.3 


< 2.7 


1 v± JJ tl ^ 


55OM 




4 


0.0 


< 33.8 


1.3 - 179.0 


< 37.7 


1.6 - 329.1 


< 39.4 


1.6 - 348.8 






0.5 


< 32.6 


< 13.0 


< 29.9 


< 7.2 


< 12.2 


< 3.8 


NGC 253 ULX2 


1 


0.0 


A f D TT > 

1 V1 ti tl <^ 


37 7Mn 


11.8 - 14.8 


413 1 - 601 


12.7 - 14.7 


471 2 - 629 8 






0.5 


< 1.7 


76.5-85.7 


< 6.8 


30.4-69.7 


M B h < 


59OM 


NGC 253 ULX3 


1 


0.0 


3.0 - 10.5 


17.8 - 290.3 


5.3 - 13.1 


17.6 - 394.3 


5.6 - 13.1 


17.6 - 399.5 






0.5 


< 15.5 


10.7 - 117.8 


< 20.5 


6 5 - 230 5 


< 36.1 


2 7 - 272 8 


NGC 253 XMM6 


1 


0.0 


< 7.2 


1.5-46.0 


< 13.6 


2.3-103.2 


< 14.2 


2.3-103.3 






0.5 


< 5.7 


< 2.3 




M B h < 


72M 






2 


0.0 


1.2 - 1.5 


22.4 - 44.6 


1.2 - 5.5 


20 1 - 9°, 5 


1.2 - 6.0 


20.1 - 100.1 






0.5 


< 13.7 


0.9-4.4 


< 6.7 


< 2.6 


M B h < 


4O5M 


M81 X-6 


1 


0.0 


1.4 - 33.5 


3 5 - 85 6 


5.1 - 40.1 


4.4 - 124.7 


5.1 - 43.0 


4.5 - 125.1 






0.5 


< 2.0 


2.1-7.8 


< 6.7* 


0.8-4.4 


< 4.4 


< 2.0 


Hoi IX X-l 


1 


0.0 


A f R TT > 

1 v± d ti ^ 


18.6M 


8.9 - 10.5 


33.3 - 63.2 


8.9 - 10.5 


33.3 - 64.2 






0.5 


< 16.0 


5.3 - 8.8 


< 27.0 


2.1 - 6.2 


< 2.3 


1.0 - 1.4 




3 


0.0 


< 13.6 


1.3 - 46.0 


< 21.9 


1.8 - 99.0 


< 24.0 


1.8 - 99.6 






0.5 


< 6.2 


< 2.5 




M B h < 


85M 




NGC 1313 X-l 


1 


0.0 


1.5 - 32.3 


3.8 - 85.5 


6.7 - 40.2 


4.8 - 31.5 


6.7-43.7 


5.0 - 31.6 






0.5 


< 29.7 


2.5-8.7 


< 26.1 


1.0-5.8 


< 3.3 


< 1.7 


NGC 1313 X-2 


1 


0.0 


1.7 - 31.0 


4.6 - 55.1 


8.0 - 40.0 


5.7 - 24.2 


8.5 - 44.1 


5.9 - 24.2 






0.5 


< 28.7 


3.3-8.4 


< 24.9 


1.3-5.6 


< 2.1 


< 1.3 


IC 342 X-l 


1 


0.0 


M D TT > 


18.9Mq 


8.3 - 12.7 


124.8 - 461.1 


8.8 - 12.7 


128.8 - 477.6 






0.5 


No constraint 


20.1-43.1 


No constraint 


7.8-25.3 


No constraint 


3.4-12.1 




2 


0.0 


1.3-39.1 


3.5-85.7 


2.5-47.5 


4.4-266.9 


2.5-51.3 


4.5-271.8 






0.5 


< 44.9 


1.8 - 12.9 


< 42.2 


< 8.1 


< 12.9 


< 3.9 


IC 342 X-6 


1 


0.0 


No constraint 


2.4-223.0 


No constraint 


2.8-259.9 


No constraint 


2.8-264.5 






0.5 


No constraint 


1.5 - 93.8 


No constraint 


< 81.8 


No rnnstraint 


< 82.2 


NGC 2403 X-l 


1 


0.0 


1.4 - 21.3 


2.9 - 85.5 


5.0 - 27.8 


3.9 - 100.6 


5.0 - 29.7 


4.0 - 100.7 






0.5 


< 19.4 


1.5-6.2 


< 14.1 


< 4.1 


< 0.8 


< 0.7 


Hoi II X-l 


1 


0.0 


> 2.6 


21.1 - 145.6 


> 10.9 


21.4 - 164.7 


> 11.3 


21 Q _ 1 68 8 






0.5 


No constraint 


15.5-39.4 


No constraint 


7.7-25.1 


No constraint 


3.2-16.0 


M83 XMM1 


1 


0.0 


No constraint 


1 5 _ 3nn o 


No constraint 


1.9 - 479.1 


Nn constraint 


1 Q _ 485 3 






0.5 


No constraint 


< 57.3 


No constraint 


< 34.1 


No constraint 


< 11.2 


NGC 5204 X-l 


1 


0.0 


9.7-61.1 


8.3-18.3 


15.5-67.6 


8.9-19.5 


17.0-71.9 


9.2-19.4 






0.5 


< 57.6 


8.0-23.0 


< 56.1 


5.3-11.6 


< 35.7 


2.2-6.9 




2 


0.0 


No constraint 


216.9-349.9 




M B h < 


2OM 








0.5 


M B h > 


15.9M 


14.7-16.1 


272.2-506.9 


14.7-16.1 


209.9-341.1 


NGC 5408 X-l 


1 


0.0 


> 2.5 


8.2-243.9 


> 7.0 


8.9-637.2 


> 7.8 


9.2-655.6 






0.5 


No constraint 


7.9-64.8 


No constraint 


4.1-36.2 


No constraint 


1.7-17.7 


NGC 3034 ULX5 


1 


0.0 


> 11.8 


> 151.7 


> 71.7 


> 209.3 


> 71.3 


> 289.7 






0.5 


< 1.6 


7.7-117.8 


< 9.3 


> 9.3 


< 9.3 


> 10.4 



Note. — a Candidate counterpart IDs, as listed in previous tables. Irradiation model output parameters as derived from fits using Galactic 
:tinction/absorption corrected magnitudes for each candidate counterpart. All stellar constraints are given in units of solar mass and radius. 
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Table 8 

la binary parameter constraints derived from intrinsic magnitudes for each 
potential ULX optical counterpart. 









±y± d ti 


1OM 


1V± ti ti 


1OOM 


1V± ti ti 


IOOOMq 


Source 


CfP lD a 


cos(7) 




















(Mn) 


(RrO 


(MrO 


{RrO 


(MrO 




Mf^r" 1 TTT Y1 
JNoL. JJ ULA1 


i 
i 


n n 
U.U 


No constraint 


O 1 1 Q1 O 
Z. i — lol.Z 


No constraint 


o.4 — zy4. / 


No constraint 


Q A OOO O 

o.4 — Z\3\3.\3 






0.5 


No constraint 


1.4-33.4 


No constraint 


< 14.8 


No constraint 


< 9.6 




2 


0.0 


< 61.4 


2.3-117.8 


< 66.7 


2.9-261.3 


< 70.3 


2.9-266.1 






u.o 


^ ft Q 1 
<v OO.l 


l.U — Z0.4 


<. ou. i 


/ 11 o 
<C 11. 


<. oy.u 


< I .Z 


NGC4190 X-l 


1 


0.0 


> 1.7 


1.4-517.4 


> 3.3 


> 2.2 


> 3.4 


> 2.3 






U.O 


No constraint 


No constraint 


No constraint 


No constraint 


No constraint 


No constraint 


NGC 253 ULX1 


1 


0.0 


No constraint 


No constraint 


No constraint 


No constraint 


No constraint 


No constraint 






0.5 


No constraint 


No constraint 


No constraint 


No constraint 


No constraint 


No constraint 




2 


0.0 


No constraint 


No constraint 


No constraint 


No constraint 


No constraint 


No constraint 






U.O 


No constraint 


No constraint 


No constraint 


No constraint 


No constraint 


No constraint 




4 


0.0 


No constraint 


No constraint 


No constraint 


No constraint 


No constraint 


No constraint 






0.5 


No constraint 


No constraint 


No constraint 


No constraint 


No constraint 


No constraint 


NGC 253 ULX2 


1 


0.0 


> 2.8 


> 11.8 


> 8.8 


> 12.3 


> 8.8 


> 12.8 






0.5 


No constraint 


> 11.2 


No constraint 


> 7.3 


No constraint 


> 3.2 


NGC 253 ULX3 


1 


0.0 


No constraint 


No constraint 


No constraint 


No constraint 


No constraint 


No constraint 






u.o 


No constraint 


No constraint 


No constraint 


No constraint 


No constraint 


No constraint 


NGC 253 XMM6 


1 


0.0 


1.2-66.0 


2.0-218.6 


1.3-71.7 


2.7-478.8 


1.3-76.0 


2.8-482.5 






U.O 


No constraint 


< zo.l 


No constraint 


^ 1 A O 

< 14. y 


No constraint 


^ Q 1 
< O.l 




2 


0.0 


No constraint 


3.0-220.2 


No constraint 


3.7-401.0 


No constraint 


3.8-420.4 






0.5 


No constraint 


1.6-52.4 


No constraint 


< 27.8 


No constraint 


< 14.0 


M81 X-6 


1 


0.0 


> 1.5 


3.5-239.6 


> 5.6 


4.4-332.7 


> 5.6 


4.5-333.7 






0.5 


< 5.8 


2.2-53.5 


< 20.2 


0.8-35.2 


< 46.5 


< 14.6 


Hoi IX X-l 


1 


0.0 


1.6-80.4 


6.2-117.8 


8.2-94.6 


7.3-187.5 


8.2-105.3 


7.8-187.5 






U.O 


<v OO.l 


A 1 OO Q 
4. / — ZZ.o 


<. oo.y 


1 Q 1/11 
l.o — 14. 1 


<. Oo.o 


<v O.O 




3 


0.0 


< 15.1 


1.3-62.2 


< 23.9 


1.9-194.8 


< 26.3 


1.9-199.6 






0.5 


< 12.6 


< 4.2 


< 5.4 


< 2.2 


M B h < 


35OM 


NGC 1313 X-l 


1 


0.0 


2.0-69.5 


5.2-84.1 


9.6-79.4 


6.2-41.1 


9.6-85.7 


6.5-41.0 






0.5 


< 66.0 


4.4-18.4 


< 64.4 


1.7-11.2 


< 33.0 


< 6.6 


NGC 1313 X-2 


1 


0.0 


> 1.7 


4.9-215.8 


> 8.2 


6.0-132.0 


> 8.3 


6.3-132.9 






0.5 


No constraint 


3.6-42.4 


No constraint 


1.4-26.4 


No constraint 


< 14.4 


IC 342 X-l 


1 


0.0 


11.3-115.1 


22.1-437.3 


12.2-114.5 


> 23.2 


12.7-112.0 


> 25.0 






0.5 


No constraint 


15.4-255.1 


No constraint 


9.7-312.5 


No constraint 


4.3-587.2 




Z 


n n 
U.U 


> 1.0 


0.1 — ZZZ.6 


> O.z 


O.U — OU ( .Z 


> o.o 


o.o — 041.0 






0.5 


No constraint 


3.9-55.5 


No constraint 


1.6-30.6 


No constraint 


< 17.1 


IC 342 X-6 


1 


0.0 


No constraint 


No constraint 


No constraint 


No constraint 


No constraint 


No constraint 






0.5 


No constraint 


No constraint 


No constraint 


No constraint 


No constraint 


No constraint 


NGC 2403 X-l 


1 


0.0 


> 1.3 


2.6-326.8 


> 5.1 


> 3.6 


> 5.2 


> 3.6 






0.5 


No constraint 


1.2-82.8 


No constraint 


< 51.5 


No constraint 


< 21.9 


Hoi II X-l 


1 


0.0 


> 2.1 


9.8-274.0 


> 10.9 


19.6-167.7 


> 11.7 


20.6-192.3 






0.5 


No constraint 


9.5-127.8 


No constraint 


3.9-86.7 


No constraint 


1.7-45.0 


M83 XMM1 


1 


0.0 


No constraint 


No constraint 


No constraint 


No constraint 


No constraint 


No constraint 






0.5 


No constraint 


No constraint 


No constraint 


No constraint 


No constraint 


No constraint 


NGC 5204 X-l 


1 


0.0 


> 2.5 


4.4-402.8 


> 8.5 


5.1-294.6 


> 8.8 


5.3-298.4 






0.5 


No constraint 


3.4-650.2 


No constraint 


> 1.3 


No constraint 


No constraint 




2 


0.0 


> 6.2 


> 12.0 


> 13.3 


> 12.4 


> 13.4 


> 12.9 






0.5 


No constraint 


> 11.7 


No constraint 


> 10.1 


No constraint 


> 4.3 


NGC 5408 X-l 


1 


0.0 


> 1.4 


> 4.3 


> 2.7 


> 5.2 


> 2.8 


> 5.3 






0.5 


No constraint 


3.3-397.2 


No constraint 


1.3-654.2 


No constraint 


No constraint 


NGC 3034 ULX5 


1 


0.0 


> 11.8 


> 151.7 


> 71.7 


> 209.3 


> 71.3 


> 289.7 






0.5 


No constraint 


> 146.9 


No constraint 


> 160.1 


No constraint 


> 57.1 
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Table 8 — Continued 







M B h 


= 10M© 


M B h 


= 1OOM 


M BH 


= IOOOMq 


Source 


C/P ID a cos(i) 


M* 
















(M ) 


(Rq) 


(M ) 


OR©) 


(M ) 


(Ho) 



Note. — As with Table [7] "Candidate counterpart ID, as listed in previous tables. All stellar constraints are given in units of solar mass 
and radius, with values derived from irradiation models fitting intrinsic X-ray and optical emission. 



Table 9 

Black hole mass, and binary parameter constraints obtained during X-ray 
irradiation model fitting (1 a constraints). 



Source 



C/P ID a Fitting 6 cos(z) 



Black hole 

mass 
constraints 



Constraints obtained for 



companion star 











(M ) 


(M ) 


CR©) 


NGC 55 ULX1 


1 


G 


0.5 


< 700 


< 13.9 


0.9- 


5.1 




2 


G 


0.5 


< 840 


< 15.1 


1.4- 


5.4 


NGC 253 ULX1 


3 


G 


0.5 


< 550 


< 12.2 


< 4 


.6 


NGC 253 ULX2 


1 


G 


0.0 


> 37.7 


11.8-14.8 


413.1- 


629.8 






G 


0.5 


< 590 


< 6.8 


30.4- 


85.7 


NGC 253 XMM6 


1 


G 


0.5 


< 72 


< 5.7 


< 2 


.3 




2 


G 


0.5 


< 405 


< 13.7 


< 4 


.4 


Hoi IX X-l 


1 


G 


0.0 


> 18.6 


8.9-10.5 


33.3- 


64.2 




3 


G 


0.5 


< 85 


< 6.2 


< 2 


.5 




3 


I 


0.5 


< 350 


< 12.6 


< 4 


.2 


IC 342 X-l 


1 


G 


0.0 


> 18.9 


8.3-12.7 


124.8- 


477.6 


NGC 5204 X-l 


2 


G 


0.0 


< 20 


No constraint 


216.9- 


349.9 






G 


0.5 


> 15.9 


14.7-16.1 


209.9- 


506.9 



Note. — Notes: "Candidate counterpart ID, as listed in previous tables. b The type of magnitudes used in fitting X-ray irradiation models: 
G - Galactic extinction/absorption corrected values used, I intrinsic values used for fitting. All stellar constraints are given in units of solar 
mass and radius, with values derived from irradiation models fitting intrinsic X-ray and optical emission. 
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